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1.0  ABSTRACT 

'  A^-prel  iminary  investigation  of  an  aeroelastically-induced  constrained  warping 
phenomenon  for  a  composite  (supermaneuver able  type)  aircraft  wing^-^as  been 
carried  out  in  this  report /-,7)ie  wing  is  analytically  modelled  as  a  straight 

L  A  *5 ' 

flat  laminated  plater  Various  forms  of  highly  simplified  aerodynamic  loads^re^ 
employed  in  the  analysis.  The  free  vibrations  and  stability  aspects^of  this 
phenomenon  are  examined  to  obtain  ^bmtf^physi cal  insights  and  to  determine 
importance  and/or  design  implications.  ^Analytical  tools  employed  include  an 
affine  transformation  concept ^vhich  was  formulated  previously  (by  the  Principal 
Investigator)  as  well  as^  non-dimensional ization  scheme C.  With  the  help  of  - 
c;  the$e~tee+&,  an  evolution  of  effective  warping  parameters  with  which  to  study 

Ut  M  ,-ri. 

this  phenomenon  was/^arrfed  out/  The  virtual  work  theorem  and  variational  prin¬ 
ciples  were  used  to  derive  the  equations  of  motion  based  on  the  assumed  wing 
displacements.  Closed  form  solutions- which  are  generally  expected  to  be 
extractable) *to  the  uncoupled  versions. of  these  equations  were  ■carefully 

A  I  *  '■>  v 

examined ^and  the  following  results l^iwre  revealed: — (d J> incorrect  modelling  of 
the  warping  phenomenon  can  lead  to  errors  in  excess  of  80%  for  the  analytically 
predicted  aeroelastic  characteristics  of  composite  aircraft  wings’, 

MiV  accurate  modelling  of  the  warping  phenomenon  is  particularly  important  for 
.(a)  wings  with  mass  coupling  or  elastic  coupling  (e.g.,  wings  aeroelastically 

(’  c  I 

tailored  using  elastic  coupling)  and higher  vibration  modesj^ii)  neglect 
of  the  warping  constraint  can  ^a^result  in  either  under-  predicted  or  ~ 

,  overpredicted  analytical  results  and  (b)'lead  to  incorrect  identification  of 
aeroelastic  divergence  modes.  The  existence  of  closed-form  free  vibrations  so¬ 
lutions  for  composite  wings  with  elastic  coupling  and  constraint  of  warping  was 
established.  ,  These  closed-form  solutions  revealed  that,  (i)  elastic  coupling 
lowers  the  coupled  frequencies  Hn  fact  a  significant  amount  of  coupling  could 


OYIBO: 8/8/86 
VT3/ENG : 17 


reduce  the  first  frequency  to  almost  zero),  (ii)  it  is  possible  to  determine 
the  origin  of  a  coupled  mode,  ( i i i )  the  mechanism  by  which  modal  transforma¬ 
tions  take  place  may  be  linked  to  some  conservative  inter-modal  energy  trans¬ 
fer.  It  has  also  been  found  that  there  is  an  interdependency  between  certain 
design  variables,  which  may  be  ordinarily  taken  as  independent,  in  the 
tailoring  technology.  Influence  of  chord-wise  bending  was  found  to  be  more  im¬ 
portant  for  bending  modes  than  twisting  modes  for  orthotropic  wings.  A  prelimi¬ 
nary  analysis  seems  to  suggest  that  the  constrained  warping  phenomenon  may  be 
captured  numerically  by  using  an  "extended  beam  element". 
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2.0  NOMENCLATURE 


co»co 

°1J 

e 


El,  GJ 
h 

(h0,a  o) 


K,  So 
kij 

L0,M0 

m0 

(dp,4p0) 


t 

u,u0 


'x,y,zj,  f x0,y0,z0) 


chordwise  integrals 

affine  space  half-chord  and  chord,  respectively 
elastic  constants 

parameter  that  measures  the  location  of  the  reference  axis 
relative  to  mid-chord 

bending  and  torsional  stiffness,  respectively 
wing  box  depth 

affine  space  bending  and  torsional  displacement, 
respectively 

elastic  coupling  and  warping  stiffness,  respectively 
elemental  stiffness  parameter 
Strouhal  number 

affine  space  running  aerodynamic  lift  and  moments, 
respectively 

affine  space  half-span  for  the  wing 
affine  space  mass  per  unit  span 

differential  aerodynamic  pressure  distributions  in  physi¬ 
cal  and  affine  space,  respectively 
time 

virtual  work  expressions  in  physical  and  affine  space, 

respectively 

flutter  speed 

physical  and  affine  space  coordinates,  respectively 
displacement  shape  functions 
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r,Li ,l2,D*,D*0, 


P,P. 


6 

w 

u 


=  generic  nondimensional ized  stiffness  parameters 
=  affine  mass  ratio  parameter 

=  Poisson  ratios  and  generalized  Poisson's  ratio,  respectively 
=  divergence  parameter 

=  affine  space  material  and  air  density,  respectively 
=  twisting  displacements 
=  displacement 
*  vibration  frequency 
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3.0  INTRODUCTION 

In  the  late  1960's  designers  began  to  investigate  the  possibility  of  exploiting 
the  directional  properties  of  composite  materials  to  improve  aeroelastic  charac 
teristics  of  lifting  surfaces.  For  example,  superfighter  designers  may 
aeroelastically  tailor  a  wing  so  that  it  deforms  to  the  optimum  camber  under 
maneuvering  loads. 

For  instance,  a  wing  can  be  tailored  so  that  its  leading  edge  will  twist  down¬ 
ward  under  the  stress  of  a  tight  turn,  thereby  decreasing  the  wing's  angle  of 
attack  and  hence  reduce  drag. 

Exploitation  of  the  directional  properties  of  composite  materials  to  solve 
aeroelastic  instability  problems  (known  as  aeroelastic  tailoring)  may  enhance 
the  performance  of  future  high  performance  military  aircraft,  particularly  the 
super-maneuverable  concepts.  This  stems  from  the  need  for  future  military 
weapon  systems  to  exhibit  high  performance  and  minimum  vulnerability  (e.g.;  min 
imum  radar  cross-sectional  areal.  These  requirements  may  lead  to  unorthodox 
aircraft  configurations  which  in  turn  result  in  unorthodox  aeroelastic  instabil 
ity  problems  (e.g.  the  X-29  primary  mode  of  instability,  called  the  'body 
freedom'  flutter). 

Perhaps  the  aeroelastic  tailoring  concept  would  have  been  discovered  much 
earlier  if  the  physics  of  anisotropic  aeroelasticity  were  more  apparent.  This 
inherent  physical  intractability  is  largely  due  to  the  existence  of  numerous 
variables,  e.g.;  flight  parameters,  several  composite  directional  properties, 
fiber  orientation  angles,  etc.,  which  are  not  even  truly  independent,  in  the 
anisotropic  system.  This  state  of  affairs  is  therefore  analogous  to  that  which 
existed  in  the  basic  rigid  body  aerodynamics  before  the  advent  of  the  similar¬ 
ity  rule  theory.  This  theory  clearly  revealed  that,  the  utilization  of 
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non-dimensional i zed  groups  of  variables,  e.g.;  Reynolds,  Mach,  Strouhal ,  Froude 
numbers  etc.,  provide  significantly  superior  physical  insight  to  the  problem 
than  the  individual  physical  variables.  The  new  methodology  that  is  being  used 
as  the  basic  tool  in  this  research  program  is  basically  the  aeroelastic  equiva¬ 
lent  of  the  aerodynamic  similarity  rule.  The  expected  superiority  of  this  new 
approach  over  the  state-of-the-art  (SOA)  counterpart,  which  utilizes  individual 
physical  variables,  especially  in  terms  of  physical  insights,  has  been 
demonstrated  in  references  1-10.  For  example,  a  high-aspect-ratio  composite 
wing  could  behave  aeroelastically  like  a  low  aspect  ratio  wing  and  vice-versa. 
Similarity  parameters  can  expose  conditions  for  which  this  might  happen.  This 
is  significant,  (for  instance)  in  the  light  of  the  important  role  played  by  the 
wing  aspect  ratio  in  the  aerodynamics  approximations  for  an  aeroelastic  anal¬ 
ysis.  This  result  may  therefore  be  suggesting  some  new  form  of  coupling  be¬ 
tween  the  elastic  and  aerodynamic  equations  in  composite  wing  aeroelastici ty. 

A  rundamental  aspect  of  this  new  methodology  has  been  used  in  studies  at  Purdue 
Uni  vers ityll  and  MIT12. 

In  this  rese  ~n,  an  investigation  of  a  wing's  spanwise  sectional  distortion 
(warping)  resulting  from  aerodynamic  forces  and  its  influence  on  the  wing's 
free  vibrations,  as  well  as  (aeroelastici  stability,  has  been  carried  out.  In 
particular,  St.  Venant's  torsional/twisting  theory,  which  is  currently  widely 
used  for  estimating  the  wing's  twisting  displacements  has  been  examined  with 
the  help  of  the  new  methodology,  to  determine  its  limitations,  when  applied  to 
wings  fabricated  of  composite  materials.  The  relevance  and  significance  of 
this  study  for  the  newly  emerging  supermaneuverable  type  aircraft  may  be  seen 
in  the  light  of  the  fact  that  (a)  supermaneuverability  is  characterized  by  high 
angle  of  attack  which  implies  high  twisting  aerodynamic  forces,  and  (b)  most 
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aircraft  designers  believe  future  aircraft  will  be  fabricated  of  40-70% 
composite  materials. 

When  the  St.  Venant's  torsional  theory  is  used  to  estimate  a  wing's  twisting 
displacement  and/or  forces,  the  fact  that  the  wing's  root  section's  distortion 
is  relatively  small  compared  to  that  of  other  sections  is  ignored.  However, 
previous  investigators  have  determined  that  such  an  unrealistic  assumption  may 
lead  to  only  little  errors  if  the  aspect  ratio  of  the  wing  is  very  high. 


% 

i 


The  research  has  shown  in  this  report  that  the  conclusions  reached  by  previous 
investigators  are  basically  true  for  wings  fabricated  of  metals  or  isotropic 
materials.  A  set  of  new  theories  is  therefore  postulated  as  shown  in  this  re¬ 
port  for  accurately  estimating  the  twisting  displacements,  vibrational 
frequencies  and  instability  boundaries  for  wings  fabricated  of  composite 
materials. 
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This  report  is  divided  into  8  sections:  sections  1,  2  and  3  show  the  abstract, 
nomenclature,  and  a  general  introduction,  respectively,  while  sections  4  and  5 
treat  the  research  objectives  and  the  status  of  the  research  effort,  respec¬ 
tively.  Section  5  outlines  the  main  findings  of  the  research  as  contributed  by 
Or.  G.A.  Oyibo,  of  Fairchild  Republic  who  was  the  Principal  Investigator  and 
Professor  T.A.  Weisshaar  of  Purdue  University  who  participated  in  the  research 
in  order  to  broaden  the  research  scope.  The  evolution  of  an  effective  warping 
parameter  with  which  to  study  the  constrained  warping  phenomena  for  wings 
fabricated  of  composite  materials  as  well  as  a  preliminary  analysis  were 
formulated  by  Dr.  G.A.  Oyibo  and  presented  in  sections  5.2  and  5.3.  These  sec¬ 
tions  basically  identify  the  important  paramters  that  physically  define  or  af¬ 
fect  these  phenomena.  A  major  accomplishment  of  this  research  program,  also 
developed  by  Dr.  G.A.  Oyibo,  was  the  formulation  of  the  closed-form  solutions, 
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for  the  first  time,  to  the  elastically  coupled  free  vibration  problems  of 
composite  aircraft  wing  undergoing  constrained  warping  and  presented  in  section 
5.4.  Professor  T.A.  Weisshaar  has  presented  the  effect  of  constrained  warping 
on  the  divergence  of  orthotropic  type  wings  in  section  5.5,  while  in  section 
5.6  he  has  presented  the  results  of  his  study  on  the  interdependence  between 
the  parameters  involved  in  aeroelastic  tailoring.  Professor  T.A.  Weisshaar  has 
also  presented  the  results  of  his  preliminary  investigation  of  the  possibility 
of  developing  a  numerical  extended  beam  superelement  that  is  capable  of 
capturing  the  essence  of  constrained  warping  in  aircraft  wings  fabricated  of 
composite  materials.  The  accomplishment  of  this  task  shall  enable  one  to  ana¬ 
lyze  wings  of  more  arbitrary  geometric  configuration. 

Section  6  provides  the  brief  overall  conclusions  and  recommendations  for  future 
work  for  furthering  the  understanding  of  these  phenomena. 


In  Sections  7,  8,  9  and  10  are  shown  status  of  publication,  the  list  of  person¬ 
nel,  interactions  with  the  Air  Force  agencies  and  discoveries  and  inventions, 
respectively. 
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4.0  RESEARCH  OBJECTIVES 

The  overall  goal  of  this  research  program  is  twofold:  first,  is  to  evolve  an 
effective  warping  parameter  from  an  advanced  lamination  theory,  with  which  to 
study  the  underlying  physics  of  the  mechanism  of  flutter  and  divergence  in 
composite  aircraft  wings  and  second,  to  use  cantilevered  flat  plate  models  of 
aircraft  wings  along  with  two-  and  three-dimensional  unsteady  aerodynamic  theo¬ 
ries  to  conduct  a  study  of  aeroelastic  stability  boundaries  and  their  depen¬ 
dence  on  the  warping  parameter. 

The  results  of  the  efforts  towards  achieving  this  overall  objective  are  the 
subject  of  this  final  report.  The  initial  three  month-effort  was  aimed  at 
conducting  investigations  to  examine  various  ways  the  warping  parameter  may  be 
evolved  and  to  determine  which  of  these  evolutions  could  best  be  utilized  to  ef¬ 
fectively  study  the  effects  of  the  "restraint  of  warping"  on  the  mechanism  of 
vibration,  aeroelastic  divergence  and  aeroelastic  flutter  of  aircraft  wings 
fabricated  of  composite  materials.  The  two  main  tools  employed  in  conducting 
these  research  studies  were  the  new  mathematical  modelling  concept  recently 
developed  by  the  Principal  Investigator,  as  well  as  a  non-dimensionalization 
scheme  which,  among  other  things,  served  as  a  tool  for  comparison.  The  major 
goal  of  the  second  quarter  was  to  research  the  possibility  of  obtaining 
closed-form  solutions  to  the  free  vibrations  problems  of  composite  aircraft 
wings  having  a  warping  constraint  and  elastic  coupling,  a  task  that  apparently 
had  never  been  successfully  accomplished  previously.  The  accomplishment  of 
this  task  is  generally  believed  to  be  vital  to  the  understanding  of  the  dynamic 
or  stability  phenomena  of  such  aircraft  wings.  Also  included  in  the  objectives 
for  the  second  quarter  were  (a)  a  preliminary  investigation  of  the  effects  of 
chordwise  bending  on  the  vibrations  of  orthotropic  plates;  (b)  examination  of 
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interdependency  of  design  variables  in  tailoring  technology,  and  fc)  the  formu-  | 

lation  of  explanations  for  some  startling  new  results  obtained  during  the  first  | 

t 

quarter.  ' 

! 

The  goals  of  the  third  and  fourth  quarters  were  (a)  to  refine  the  newly  i 

discovered  closed-form  solution  approach  in  order  to  use  it  to  study  the  ! 

( 

constrained  warping  phenomena  at  higher  vibrational  frequencies,  (b)  to  assess  J 

I 

the  implications  and  impact  of  what  has  been  accomplished  so  far  on  the 
direction  of  this  research  program,  and,  (c)  to  prepare  the  AIAA  Paper  No.  86- 
1006,  a  paper  based  on  the  preliminary  results  of  this  research  program,  which 
was  presented  at  the  AIAA/ASME/ASCE/AHS  27th  Structures,  Structural  Dynamics 
and  Materials  Conference,  held  in  San  Antonio,  Texas,  May  19-21,  1986. 
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I 

5.0  STATUS  OF  RESEARCH  EFFORT  j 

I 

5.1  INTRODUCTION  ! 

t 

During  the  first  year,  the  research  program  progressed  according  to  plan  j 

l 

resulting  in  the  accompl ishment  of  the  goals  defined  for  the  period.  These 
goals,  which  have  been  defined  in  Section  2  of  this  report,  basically  include 
a  preliminary  formulation  and  investigation  of  the  aeroelasticity  of  supermaneu¬ 
ver  able  type  aircraft  wings  (fabricated  of  composite  materials)  in  order  to  de¬ 
termine  an  effective  and  efficient  means  with  which  to  accurately  estimate  and 
study  the  effects  of  the  "restraint  of  warping"  phenomenon  on  their  vibration 
and  mechanism  of  aeroelastic  instability. 

In  the  first  quarter,  during  the  preliminary  study.  Dr.  6. A.  Oyibo  investigated 
the  free  vibrations  and  the  aeroelastic  flutter  aspect  of  this  phenomenon, 
using  an  affine  transformation  concept  and  Professor  T.A.  Weisshaar  studied  the 

aeroelastic  divergence  aspect  (of  the  phenomenon),  using  a  non- 
dimensionalization  scheme.  The  results  obtained  seem  to  indicate  that  although 
three  different  problems  were  investigated  (using  two  different  approaches) 
there  were  some  general  phenomenological  trends  that  appeared  to  be  common  to 
all  three  of  them.  Summarized  versions  of  these  investigations  were  presented 
in  the  first  quarterly  report.  Journal  articles  based  on  the  results  of  these 
investigations  are  being  prepared  to  be  submitted  to  the  AIAA  Journal  of  Air¬ 
craft  for  publication. 

In  the  beginning  of  the  second  quarter.  Dr.  Oyibo  was  able  to  formulate  an  ex¬ 
planation  for  the  startling  new  results  obtained  during  the  first  quarter.  Dur¬ 
ing  the  second  quarter  Dr.  Oyibo  was  also  able  to  establish  the  existence  of 
closed-form  solutions  for  the  free  vibrations  of  composite  aircraft  wings  with 
the  constraint  of  warping  and  elastic  coupling.  The  preliminary  results,  which 
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appear  to  be  exciting,  confirm  some  known  trends  and  expose  some  new  ones.  Pro¬ 
fessor  T.A.  Weisshaar  examined  the  interdependence  between  some  design  vari¬ 
ables  employed  in  the  tailoring  technology  as  well  as  the  effects  of  chordwise 
bending  on  the  frequencies  of  orthotropic  plates. 

A  more  in-depth  investigation  of  the  newly  discovered  closed-form  solutions  was 
carried  out  during  the  third  and  fourth  quarters  to  further  investigate  the  dy¬ 
namic  phenomena  of  composite  aircraft  wings  with  the  constraint  of  warping  and 
elastic  coupling. 

In  particular,  higher  vibrational  frequencies  are  being  generated  and  studied 
carefully  to  determine  trends.  An  assessment  of  some  of  the  impacts  and  impli¬ 
cations  of  results  obtained  in  the  research  thus  far  reveals  that  (a)  the  newly 
discovered  closed-form  approach  may  be  extendable  to  study  the  stability  prob¬ 
lems  and  (b)  it  now  appears  possible  (and  certainly  desirable)  to  develop  a  nu¬ 
merical  superelement  that  can  capture  the  essence  of  the  warping  phenomena  for 
studying  more  irregular  wing  configurations.  During  the  third  quarter,  also, 
AIAA  Paper  No.  86-1006,  based  on  the  preliminary  results  of  this  research  pro¬ 
gram,  and  accepted  to  be  presented  at  the  AIAA/ASME/ASCE/AHS  27th  Structures, 
Structual  Dynamics  and  Materials  Conference  in  San  Antonio,  Texas,  May  19-21, 
1986  was  prepared,  and  subsequently  presented  at  the  conference. 
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5.2  ACCURATE  DYNAMIC  THEORY  FOR  SUPERMANEUVERABLE  AIRCRAFT  WINGS 
by  Gabriel  A.  Oyibo 

5.2.1  Introduction 

Modern  supermaneuverable  aircraft  concepts  benefit  a  great  deal  from,  among 
other  things,  significant  advances  in  materials  technology  and  the  availability 
of  more  accurate  aerodynamic  prediction  capabilities.  Supermaneuverability  as 
a  design  goal  invariably  calls  for  an  optimization  of  the  design  parameters. 
Optimization  may  be  partially  accomplished  for  example,  by  using  composite  mate¬ 
rials  to  minimize  weight.  Indeed,  it  has  been  known  that  these  composite  mate¬ 
rials  can  be  tailored  properly  to  resolve  the  dynamic  or  static  instability 
problems  of  these  types  of  aircraft.  The  concept  is  referred  to  as  aeroelastic 
tailoring. 

While  aeroelastic  tailoring  has  tremendous  advantages  in  the  design  of  an  air¬ 
craft,  the  analysis  which  provides  the  basis  for  the  aeroelastic  tailoring  it¬ 
self  is  generally  very  involved.  This  is  rather  unfortunate  since  a  good  funda¬ 
mental  physical  insight  of  the  tailoring  mechanism  is  required  for  accurate  and 
reliable  results. 

In  this  investigation  an  attempt  is  made  to  look  at  some  dynamics  theories  that 
can  be  used  to  understand  the  aeroelastic  tailoring  mechanism.  Specifically, 
the  accuracy  of  the  St.  Venant  torsion  theory  which  is  relatively  simple  and 
frequently  used  in  aeroelastic  analysis  is  examined  with  particular  reference 
to  the  effects  of  the  wings  aspect  ratio  as  well  as  other  design  parameters. 

An  accurate  torsion/twist  theory  is  particularly  significant  for  super¬ 
maneuverable  aircraft  wings  since  supermaneuverability  is  basically 
characterized  by  high  angle  of  attack. 
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Although  earlier  studies(l,2,3)  have  indicated  that  the  St.  Venant's  torsion 
theory  is  reasonably  accurate  except  for  aircraft  wings  with  fairly  low  aspect 
ratios,  the  theory  supporting  that  conclusion  was  based  on  the  assumption  that 
the  wing  is  constructed  of  isotropic  materials.  Basically  the  St.  Venant's  tor 
sion  theory  assumes  that  the  rate  of  change  of  the  wing's  twist  angle  with  re¬ 
spect  to  the  spanwise  axis  is  constant.  This  assumption  is  hardly  accurate  par 
ticularly  for  modern  aircraft  construction  in  which  different  construction  mate 
rials  are  employed  and  the  aerodynamic  loads  vary  significantly  along  the 
wing's  span.  However,  References  1-3  have  shown  that  (in  spite  of  such  an 
inaccurate  assumption)  the  main  parameter  that  determines  the  accuracy  of  the 
St.  Venant's  theory  is  the  wing's  aspect  *,atio.  Thus,  it  was  determined  that 
the  theory  is  fairly  accurate  for  moderate  to  high  aspect  ratio  wings  con¬ 
structed  of  isotropic  materials.  In  recent  studies  (4,5,6)  however,  it  has 
been  shown  that  for  wings  constructed  of  orthotropic  composite  materials,  the 
conclusion  of  References  1-3  need  to  be  modified.  Rather  than  using  the  geomet 
ric  aspect  ratio  of  the  wing  to  determine  the  accuracy  of  St.  Venant's  twist 
theory,  it  was  suggested  that  a  generic  stiffness  ratio  as  well  as  an  effective 
aspect  ratio  which  considers  the  wing's  geometry  and  the  ratio  of  the  principal 
directional  stiffness  should  be  considered  in  establishing  the  accuracy  of  St. 
Venant's  theory. 

The  present  investigation  is  related  to  the  studies  that  were  initiated  in  Ref¬ 
erences  5  and  6.  In  this  study  the  first  task  was  to  examine  the  role  of 
coupling  (both  mass  and  elastic  coupling)  on  the  accuracy  of  St.  Venant's  the¬ 
ory  applied  to  static  problems.  It  was  discovered  that  coupling  plays  a  very 
significant  role  on  the  accuracy  of  St.  Venant's  twist  theory.  The  second  task 
was  to  investigate  the  torsional  vibration  for  a  flat  plate  model  of  an  air¬ 
craft  wing  fabricated  of  composite  materials  in  which  the  constrained  warping 
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phenomenon  is  more  realistically  represented,  with  particular  emphasis  on 
higher  frequencies  and  to  compare  results  with  those  from  a  representation 
based  on  St.  Venant's  theory. 

5.2.2  Formulation 

Consider  an  aircraft  wing  fabricated  of  composite  materials  and  mathematically 
idealized  as  a  cantilevered  plate  subjected  to  an  aerodynamic  flow  over  its  sur 
faces.  The  mathematical  statement  of  the  virtual  work  theorem  for  such  a  plate 
model  is  well  known  and  documented.  It  is  also  known  that  such  mathematical 
statements  of  the  virtual  work  theorem  for  a  laminated  plate  model  are  char¬ 
acterized  with  the  existence  of  so  many  variables  (in  the  statement), 
reflecting  the  various  directional  properties  for  the  laminated  plate  model, 
which  would  tend  to  interfere  with  any  physical  insight  that  might  be  desired 
from  a  phenomenological  analysis  employing  such  a  mathematical  statement.  The 
newly  discovered  affine  transformation  concept  (5,6  and  7)  was  developed 
principally  to  resolve  such  a  problem. 

This  new  concept  therefore  can  be  used  to  evolve  the  mathematical  statement  of 
the  virtual  work  theorem  in  an  affine  space  given  by  the  following  equation. 
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Dij  are  the  elastic  constants,  P  ,  is  the  material  density,  dp  is  the  differ¬ 
ential  pressure  distribution,  w  is  the  displacement,  t  is  the  time,  A  integrals 
represent  area  integrals  and  h  is  the  wing  box  depth. 

Equations  (1)  and  (2)  therefore  form  the  basis  of  the  newly  developed  method¬ 
ology.  The  equations  of  motion  of  a  plate  model  of  an  aircraft  wing  can  now  be 
derived  by  prescribing  a  realistic  wing  displacement  and  using  Equation  (1). 

Xhen  Equation  (1)  is  compared  to  its  physical  space  counterpart,  it  is  seen 
that  Equation  (1)  has  fewer  variables.  It  is  also  seen  that  Equation  (1)  con¬ 
tains  only  non-dimensional ized  stiffness  quantities  (compared  to  dimensional 
stiffness  quantities  in  its  physical  space  counterpart) .  Another  feature  of 
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this  new  methodology  which  makes  it  unique  is  that  the  non-dimensional ization 
(a  consequence  of  the  affine  transformation)  is  accomplished  before  assuming 
the  wing  deformations.  This  means  that  the  non-dimensional ization  is  independ¬ 
ent  of  how  the  wing  deforms.  A  non-dimensional ization  scheme  that  depends  on 
a  particular  assumption  of  the  wings  deformations  could  lead  the  analyst  to  an 
incorrect  physical  interpretation  of  results,  since  the  wing's  deformations  as¬ 
sumptions  have  inherent  errors  because  they  are  based  on  the  analyst's  judg¬ 
ments  and  experience.  This  observation  may  become  clearer  during  the  evolution 
of  a  warping  parameter  with  which  to  study  the  effects  of  the  warping  con¬ 
straint  phenomenon  on  the  status  and  dynamics  of  a  wing  fabricated  of  composite 
materials  later  in  this  study. 

If  the  chordwise  curvature  is  neglected  in  an  initial  approximation,  the  wing's 
deformation  may  be  assumed  as  follows: 

w(t,xQ,yo)  =  h0(t.y0)  +  x0ao^t,yo^  ^ 

where  h0  and  a  0  are  the  bending  and  twisting  displacements  respectively. 

It  can  be  shown  that  when  Equation  (3)  is  substituted  into  Equation  (1)  and  the 
variational  calculus  is  carried  out  for  arbitrary  h0  and  a  0,  the  following 
equations  of  motion  are  obtained: 
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5.2.3  Evolution  of  Warping  Parameters 


The  evolution  of  the  warping  parameter  with  which  to  study  the  aeroelastic  warp 


ing  constraint  phenomenon  for  wings  fabricated  of  composite  materials  is  a  pro¬ 


cess  that  depends  on  the  sophistication  of  the  wing's  mathematical  model; 


whether  coupling  effects  are  included,  whether  the  wing's  chordwise  curvatures 


are  included  and  so  on.  Therefore,  any  warping  parameter  is  as  good  as  the  cot" 


responding  wing's  displacements  assumptions.  However,  Equation  (1)  makes  it 


possible  for  the  analyst  to  determine  its  effective  independent  variables  even 


before  the  displacement  assumptions  are  made. 


K  t* 
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By  non-dimensionalizing  the  spanwise  space  variable  in  Equation  (4),  depending 
on  whether  one  is  interested  in  the  static,  dynamic,  coupled  or  uncoupled  dis¬ 
placements,  one  of  the  following  warping  parameters  may  be  useful. 


Xc  =  ^ 


(7) 
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(/0/c0)  is  defined  as  the  wing's  effective  aspect  ratio  and  D0*  and  L  are  the 
generalized  stiffness  and  coupling  ratios  respectively  (defined  in  earlier 
work  such  as  References  5  and  6). 

Equations  (7)  thru  (9)  represent  the  appropriate  warping  parameter  for  dynamic 
deformation,  static  displacement  with  elastic  cross-coupling,  and  static  defor¬ 
mation  with  "geometric"  coupling  (  e  i6  1)  respectively. 

It  was  discovered  in  this  study  that  evolving  the  warping  parameter  in  a  manner 
shown  in  Equations  (1)  thru  (3),  should  enable  one  to  investigate  the  effects 
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of  warping  on  the  composite  wing's  dynamics  (or  the  accuracy  of  St.  Venant's 
theory!  effectively.  From  the  lamination  theory  for  composites  it  is  known 
that  while  0o*  and  f^0/c0)  are  always  positive,  L  and  a?  can  &e  positive  or  neg¬ 
ative.  However,  from  Equations  (1)  and  (3),  it  is  clear  that  whether  a  com¬ 
posite  wing  has  positive  or  negative  coupling,  the  warping  effect  (in  terms  of 
X c  )  is  unchanged. 

5.2.4  Results  and  Conclusions 


By  using  the  evolved  warping  parameters  defined  in  Equations  (7)  and  (8)  and 
the  appropriate  boundary  conditions,  the  boundary  value  or  the  eigenvalue  prob¬ 
lems  associated  with  Equation  4  are  solved  in  a  closed-form  manner  to  determine 
the  wing's  static  twist  and  torsional  vibrations  respectively. 


The  results  are  shown  in  Figures  1  thru  10.  Figure  1  shows  a  comparison  of 
the  static  wing  tip  twist  obtained  in  the  presefit  study  and  that  obtained  via 
St.  Venant's  twist  theory  in  the  presence  of  statically  distributed  aerodynamic 
forces  and  low  to  moderate  coupling.  Figures  2  thru  5  show  the  trend  for  con¬ 
centrated  aerodynamic  forces  and  substantial  coupling.  In  Figures  1  thru  5 
it  is  seen  that  the  presence  of  coupling  makes  the  errors  of  St.  Venant's  the¬ 
ory  worse.  This  seems  to  suggest  that  the  more  sophisticated  theory  is  more 
important  for  wings  with  coupling  (e.g.,  wings  aeroelastically-tai lored  using 
elastic  cross-coupling). 


Figures  2  thru  5  also  show  that  nonconservative  errors  (  )a  m. 

STV  * 

In  Figures  6  thru  10  the  results  of  the  torsional  vibrations  are  shown  for  a  wing 


(I'wi’H 


>1)  are  possible 


without  coupling.  Figure  6  is  particularly  striking  by  showing  clearly  that 
the  errors  of  St.  Venant's  theory  are  more  significant  for  higher  modes  of 
torsional  vibrations.  While  the  errors  could  be  bad  enough  for  the  first  mode 
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for  Ac  <4,  it  is  encouraging  to  find  that  errors  of  the  lower  modes  are  better 
than  those  for  higher  modes  since  the  critical  aeroelastic  instabilities  generally 
occur  at  lower  frequencies.  In  Figures  7  thru  10  the  frequency  parameters 
are  plotted  against  the  effective  aspect  ratio  of  the  wing  for  each  Oq  for 
the  first  four  modes.  The  St.  Venant's  results  should  be  recovered  for  very 
large  effective  aspect  ratio. 

Using  Figures  1  thru  10,  the  following  conclusions  can  be  summarized:  (i) 
ignoring  warping  arbitrarily  using  St.  Venant's  theory  could  result  in  very  sig¬ 
nificant  errors  (as  high  as  over  80%  errors)  in  analytical  results  for 
composite  aircraft  wings,  (ii)  warping  is  more  important  (St.  Venant's  theory 
is  less  accurate)  for  wings  with  coupling,  (iii)  St.  Venant's  theory  (which  has 
always  been  shown  to  be  conservative  (1,2)),  can  be  non-conservative  or  St. 
Venant's  approximation  can  lead  to  an  unsafe  design  error  (under  design  rather 
than  over  design  from  a  stability  point  of  view),  and  (iv)  warping  is  more  im¬ 
portant  (St.  Venant's  theory  is  less  accurate)  for  higher  vibration  modes. 
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5.3  AEROELASTIC  FLUTTER  PHENOMENON  FOR  COMPOSITE  AIRCRAFT  WINGS  WITH  WARPING 
CONSTRAINT  by  Gabriel  A.  Oyibo 

5.3.1  Introduction 

As  one  goes  from  aeroelastic  research  and  development  to  applications,  it  is  ab 
solutely  important  that  one  clearly  understands  the  limitation  of  the  theoreti¬ 
cal  idealizations  employed  during  the  research  and  development  stage.  Such  the 
oretical  idealizations  may  have  intricate  practical  implications,  particularly 
for  aircraft  lifting  surfaces  which  are  fabricated  of  composite  materials,  due 
to  the  inherent  mathematical  and  physical  understanding  difficulties  generally 
encountered  in  the  analysis  of  composite  systems. 

The  significance  of  the  above  mentioned  point  may  be  evaluated  in  the  following 
context:  it  has  been  established(l,2,3)  that  the  effects  of  spanwise  sectional 
warping  constraints  on  the  aeroelastic  phenomena  like  the  bending-torsion  flut¬ 
ter  may  be  ignored  for  a  warping  constrained  cantilevered  lifting  surface  pro¬ 
vided  that  tilt;  surface  has  a  high  structural  aspect  ratio.  A  "freely  warping" 
lifting  surface  in  this  context  may  be  defined  as  one  whose  twist  distribution 
varies  linearly  with  its  spanwise  length  (measured  from  its  root).  Therefore, 
since  it  is  almost  impossible  in  practice  to  construct  aircraft  wings  with  such 
a  twist  distribution,  it  is  probably  safe  to  assume  that  most  wings  have 
constrained  warping.  The  idealization  of  wings  as  "freely  warping",  which 
simplifies  aeroelastic  analysis  considerably,  is  based  on  the  conclusions  from 
investigations  like  those  of  References  1-3,  in  which  it  was  shown  that  a  high 
aspect  ratio  wing  may  be  treated  aeroel astical ly  as  "freely  warping".  However, 
as  was  pointed  out  in  References  4  and  5,  these  conclusions  are  strictly  true 
for  lifting  surfaces  constructed  of  isotropic  materials,  and  therefore  should 
be  modified  (or  generalized)  for  those  fabricated  of  composite  materials.  The 
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suggested  general izations(4, 5)  were  basically  meant  to  caution  analysts  from 
using  idealizations  which  were  strictly  meant  for  lifting  surfaces  fabricated 
of  isotropic  materials  for  those  fabricated  of  composite  materials. 


The  study  of  the  significance  of  the  warping  constraint  in  the  binary  aero- 
elastic  flutter  phenomenon  for  orthotropic  composite  wings  is  the  subject  of 
this  section.  In  Reference  6,  the  study  of  the  binary  flutter  mechanism  for 


orthotropic  composite  wings  was  carried  out  ignoring  the  warping  constraint 
(permitting  the  use  of  St.  Venant's  twist  theoryl.  The  results  of  the  present 
study  which  considers  the  warping  constraint  in  the  analysis  is  eventually 
compared  to  those  of  Reference  6  and  a  set  of  interesting  conclusions  is  drawn. 
Two-dimensional  unsteady  aerodynamics  and  an  advanced  composite  lamination  the¬ 
ory  were  used  in  both  studies  along  with  a  new  dynamic  equivalence  concept. (6) 
This  new  concept  is  based  on  the  following  newly  defined  generic  quantities. 


Consider  a  composite  wing  with  an  effective  aspect  ratio,  i*0/c0,  given  by 

(L'Xttf"? 


and  generic  stiffness  ratios  D*  and  Do*  given  by. 


D*  = 


°12  4  2D66  •  • 
-  ;  0o  -  0  (1  -«), 


(Du022) 

where  €  is  a  generalized  Poisson's  ratio  effect. 


5.3.2  Formulation 


The  traditional  approach  for  formulating  the  equations  of  motion  for  flat  plate 
models  of  aircraft  wings  fabricated  of  composite  materials  is  well  documented 
and  known  to  lead  to  mathematical  equations  containing  numerous  variables  with 
no  known  bounds.  The  resulting  flutter  analysis  is  generally  characterized 
with  results  that  are  hardly  explicable.  The  new  methodology,  which  provides 
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solutions  to  such  problems,  employs  an  affine  transformation  concept  to  define 
equivalent  dynamic  systems  with  which  to  investigate  the  mechanisms  of  flutter. 
Such  dynamically  equivalent  systems  automatically  provide  easy  means  for  de¬ 
signing  wind  tunnel  test  models  which  are  accurate.  The  new  methodology  basi¬ 
cally  transforms  the  mathematical  statement  of  the  virtual  work  theorem  for  a 
laminated  plate  model  of  an  aircraft  wing,  undergoing  aeroelastic  oscillations, 
into  a  nore  convenient  counterpart  in  an  affine  space,  based  on  Figure  1,  given 
by  the  following  equation. 
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Dij  are  the  elastic  constants,  p  ,  is  the  material  density,  Ap  is  the  dif¬ 
ferential  pressure  distribution,  w  is* the  displacement,  t  is  the  time,  A  inte¬ 
grals  represent  area  integrals  and  h  is  the  wing  box  depth. 

Equations  (3)  and  (4)  represent  the  basis  for  the  new  methodology.  By  pre¬ 
scribing  realistic  aircraft  wing  deformations,  these  equations  can  be  used  to 
derive  equations  of  motion. 

5.3.3  New  Similarity  Theory 

Since  the  equations  of  motion  for  the  aircraft  wing  are  derivable  from  the 
virtual  work  theorem  given  by  Equation  (3),  it  is  now  clear  (even  before  the 
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actual  equations  of  motion  are  derived)  that  the  parameters,  whose  similarities 
must  be  preserved  between  the  real  aircraft  wing  and  its  model  are  those  pre¬ 
sent  in  Equation  (3).  Therefore,  rather  than  trying  to  preserve  the  similari¬ 
ties  of  Djj,  as  well  as  the  fiber  orientation  angles  as  may  be  necessary  in 

Equation  (1),  Equation  (3)  makes  it  possible  to  require  only  the  preservation 

°il 

of  the  similarities  of  D0*,  (  -r1 —  ),  Lj,  and  L2  over  and  above  the  require- 

22 

ments  for  aircraft  wings  fabricated  of  isotropic  materials.  In  what  follows  a 
specific  displacement  assumption  will  be  used  to  illustrate  this  new  concept. 

Consider  the  following  bending-torsion  displacement  assumption: 

w(t.xo.y0)  ■  'Vt,yo^  +  x0ao(t*yo)  (5) 


where 


h0  and  a0  are  the  bending  and  torsional  displacements  respectively. 

When  Equation  (5)  is  substituted  into  Equation  (3)  the  following  equations  of 
motion  for  the  wing  can  be  derived  with  the  help  of  variational  calculus. 
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Equations  (5)  therefore  give  the  coupled  equations  of  motion  for  a  flat  plate 
model  of  an  aircraft  wing,  fabricated  of  composite  materials  undergoing  aero- 
elastic  oscillations.  The  presence  of  aj  in  Equation  (6)  indicates  a  realistic 
representation  of  aeroelastic  oscillations  with  constrained  warping.  St. 

Ven ant's  theory  which  ignores  a3  will  be  compared  to  the  present  theory  in  the 
flutter  analysis  that  follows. 


5.3.4  Unsteady  Aerodynamics 

In  order  to  solve  Equations  (6)  to  determine  aeroelastic  instability  bound¬ 
aries,  the  aerodynamic  forces  and  moments  L0  and  M0  would  have  to  be  generated 
by  solving  the  aerodynamic  equations.  In  the  affine  space,  the  aerodynamic 
equations  are  slightly  more  difficult  to  solve  than  in  the  physical  space.  How¬ 
ever,  the  overall  relative  simplicity  of  the  aeroelastic  equations  of  motion 
certainly  justifies  the  extra  effort  needed.  Two-dimensional  linear  approxima¬ 
tions  and  a  Theodorsen  type  approach  can  be  used  to  generate  the  incompressible 
harmonic  unsteady  aerodynamics  used  in  this  analysis.  It  can  be  shown  that: 


L  *  —  ^  c/  lh  h  +  l  -  L.  (e  *  y) 
o  r  oo  o  h0  0  °0  h0  1 
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L  ^ 


<Loo‘T>(e4  ?)]<*. 


40 


OYIBO: 2 :8/8/86 
VT3/ENG:17 


- 

»o  J2  2 


°22  1/4 
r  *  (  — ) 

°11 


(9) 


tho*«o)  *  (Vao)e 


Iwt 


jO  is  the  air  density,  (H0,  c(  0)  are  the  vibration  amplitudes,  u  is  the 
vibration  frequency  and  e  measures  the  location  of  the  elastic  axis  relative 
to  the  mid-chord. 


Lh0  is  the  running  aerodynamic  lift  coefficient  due  to  the  wing's  bending  vibra¬ 
tions  in  affine  space,  La  is  the  running  aerodynamic  lift  coefficient  due  to 
the  wing's  torsional  vibrations  about  the  affine  space  quarter-chord.  MeQ  is 
the  running  aerodynamic  moment  coefficient  due  to  the  wing's  torsional  vibra¬ 
tions  about  the  affine  space  quarter-chord.  These  coefficients,  which  are  func¬ 
tions  of  the  affine  space  reduced  frequency  (Strouhal  number),  k0,  are  expres¬ 
sed  as  follows: 


l 


1  i 


[l  ♦  2f(k0)] 


M 


(10) 
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where: 


C  (k0)  is  the  Theodorsen's  function  in  the  affine  space  given  by: 

CIO 


’  K<ko)  *  (H0(ko>] 


(11) 


rf- 0  and  H^  are  Hankel  functions  and  k 0  is  given  by: 


U0,  the  affine  space  flight  velocity,  is  given  by: 


(121 


UQ  *  rU  (13) 

The  Hankel  functions  in  Equation  (11)  are  usually  available  in  tabulated  forms. 
Hence  L h0»  La0»  and  MaQ,  can  &e  calculated  for  a  chosen  value  of  k0  from 
which  l0  and  M0  could  be  evaluated  for  flutter  analysis. 

It  must  be  pointed  out  that  the  affine  space  unsteady  aerodynamics  is  not 
restricted  to  two-dimensional  space  (airfoil  section),  but  could  also  be 
evaluated  in  three-dimensional  space  with  various  levels  of  sophistication 
(including  nonlinear  treatments). 

5.3.5  Flutter  Analysis 

Without  any  loss  of  generality,  Equation  (6)  shall  be  specialized  here  to  deal 
with  the  flutter  arising  from  the  coalescence  between  the  first  bending  and  the 
first  torsional  frequencies  of  the  wing.  Hence  Equations  (6),  no  matter  what 
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order  of  approximation  was  used  for  estimating  the  uncoupled  natural  frequencies. 


reduce  to: 
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rcfQ  is  basically  a  non-dimensionalized  radius  of  gyration  in  the  affine  space. 

Using  the  simple  harmonic  motion  assumptions  in  Equation  (9),  non- 
dimensionalized  closed  form  flutter  equations  resulting  from  Equations  (8)  and 


(13)  are: 
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where: 


( .597)‘yr 

r 77 


( .597  Jr)' 


where: 


Equations  (18)  and  (19)  give  the  St.  Venant's  approximation  and  the  present 
approximation  respectively.  Clearly  Equation  (16),  which  can  be  put  in  the 
form, 


are  eigenvector  equations.  Hence  for  nontrivial  solutions  the  determinant  of 
M  must  vanish  identically.  The  fact  that  M  is  complex  (resulting  from  the 
non-dimensional ized  unsteady  aerodynamic  functions)  makes  the  computations  of 
the  eigenvalues  for  a  set  of  configuration  parameters  a  tedious  exercise. 

The  computations  were  carried  out  by  selecting  sets  of  configuration  parameters 
and  solving  for  k0  ro0,M0,  *a0  and  w  •  These  parameters  are  the  affine  space 
quantities. 
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where: 
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5.3.6  Results  and  Conclusions 

The  computations  of  this  flutter  analysis  to  study  the  effects  of  proper  repre 
sentation  of  the  wing's  constrained  warping  phenomenon  is  still  in  progress. 
However  the  results  obtained  thus  far  which  appear  to  be  interesting,  permit 
the  following  conclusions  or  observations  on  the  mechanism  of  flutter  for  flat 
plate  models  of  aircraft  wings  fabricated  of  composite  with  a  more  realistic 
(or  accurate)  representation  of  the  constrained  warping  phenomenon. 


If  Uf  and  Uf^tv  are  the  flutter  speeds  predicted  by  the  present  more  sophisti¬ 
cated  theory  and  St.  Venant  theory  (which  neglects  warping  constraint)  respec¬ 
tively,  and  , u  *nd  are  the  radius  of  gyration,  mass  and  static  mass 

0  O  0 

unbalance  ratios  respectively.  Figures  2  thru  6  can  be  used  to  draw  the  follow¬ 
ing  conclusions  for  the  range  of  parameters  investigated:  (i)  significant 
errors  (conservative  and  nonconservative)  in  flutter  speeds  are  possible  if 
warping  is  neglected  in  composite  wing  flutter  analysis,  (ii)  the  errors  de¬ 
crease  with  the  generic  design  variable  Oo*  (unless  there  is  a  low  static  mass 
unbalance)  and  an  effective  aspect  ratio  /0/c0,  ( i i i )  flutter  speed  errors  due 
to  neglect  of  warping  constraint  tend  to  be  greater  for  a  system  with  less 
static  mass  unbalance,  (iv)  while  conservative  flutter  speed  errors 
( Uf /Uf g <  1)  appear  to  be  the  trend  in  systems  with  high  static  mass 
unbalance,  nonconservative  flutter  speed  errors  (Uf/Uf<jyV  >  1)  are  possible  in 
systems  with  low  static  mass  unbalance. 
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5.4  CLOSED  FORM  DYNAMIC  SOLUTIONS  FOR  SUPERMANEUVERABLE  AIRCRAFT  WINGS 
by  Gabriel  A.  Oyibo 

5.4.1  Introduction 

The  performance  of  modern  supermaneuverable  aircraft  can  be  made  to  benefit  a 
great  deal  from  significant  advances  in  materials  technology  and  the  availabil¬ 
ity  of  more  accurate  aerodynamic  prediction  capabilities.  Supermaneuver¬ 
ability,  as  a  design  goal,  invariably  calls  for  an  optimization  of  the  design 
parameters.  Optimization  may  be  partially  accomplished,  for  example,  by  using 
composite  materials  to  minimize  weight.  Indeed,  it  is  known  that  composite 
materials  can  be  tailored  to  resolve  the  dynamic  or  static  instability  problems 
of  these  types  of  aircraft.  This  concept  is  referred  to  as  aeroelastic  tailor¬ 
ing. 

While  aeroelastic  tailoring  can  provide  tremendous  advantages  in  the  design  of 
an  aircraft,  the  analysis  which  provides  the  basis  for  the  aeroelastic  tailor¬ 
ing  itself  is  generally  very  complex.  This  is  rather  unfortunate,  since  a  good 
fundamental  physical  insight  of  the  tailoring  mg  mechanism  is  required  for  accu¬ 
rate  and  reliable  results  and  earlier  theories  do  not  provide  sufficient 
i nsight. 

Although  earlier  studies(l,2,3)  have  indicated  that  the  St.  Venant's  torsion 
theory  is  reasonably  accurate  except  for  aircraft  wings  with  fairly  low  aspect 
ratios,  the  theory  supporting  that  conclusion  was  based  on  the  assumption  that 
the  wing  is  constructed  of  isotropic  materials.  Basically,  the  St.  Venant's 
torsion  theory  assumes  that  the  rate  of  change  of  the  wing's  twist  angle  with 
respect  to  the  spanwise  axis  is  constant.  This  assumption  is  hardly  accurate 
particularly  for  modern  aircraft  construction  in  which  different  construction 
materials  are  employed  and  the  aerodynamic  loads  vary  significantly  along  the 
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wing's  span.  However,  References  1-3  have  shown  that  (in  spite  of  such  an 
inaccurate  assumption)  the  main  parameter  that  determines  the  accuracy  of  the 
St.  Venant's  theory  is  the  wing's  aspect  ratio.  Thus,  it  was  determined  that 
the  theory  is  fairly  accurate  for  moderate  to  high  aspect  ratio  wings  con¬ 
structed  of  isotropic  materials.  In  recent  studies(4,5,6)  however,  it  has  been 
shown  that  for  wings  constructed  of  orthotropic  composite  materials,  the  conclu¬ 
sion  of  References  1-3  needs  to  be  modified.  Rather  than  using  the  geometric 
aspect  ratio  of  the  wing  to  determine  the  accuracy  of  St.  Venant's  twist  the¬ 
ory,  it  was  suggested  that  a  generic  stiffness  ratio  as  well  as  an  effective 
aspect  ratio  which  considers  the  wing's  geometry  and  the  ratio  of  the  principal 
directional  stiffness  should  be  considered  in  establishing  the  accuracy  of  St. 
Venant's  theory. 

In  Reference  7,  an  attempt  was  made  to  look  at  some  dynamic  theories  that  can 
be  used  to  understand  the  aeroelastic  tailoring  mechanism.  Specifically,  the 
accuracy  of  the  St.  Venant  torsion  theory  which  is  relatively  simple  and  fre¬ 
quently  used  in  aeroelastic  analysis  is  examined  with  particular  reference  to 
the  effects  of  the  wing's  aspect  ratio,  as  well  as  other  physical  parameters. 

In  the  study  (7)  the  first  task  was  to  examine  the  role  of  coupling  (both  mass 
and  elastic  coupling)  on  the  accuracy  of  St.  Venant's  theory  applied  to  static 
problems.  It  was  discovered  that  coupling  plays  a  very  significant  role  on  the 
accuracy  of  St.  Venant's  twist  theory.  The  second  task  was  to  investigate  the 
torsional  vibration  for  a  flat  plate  model  of  an  aircraft  wing  fabricated  of 
composite  materials  in  which  the  warping  phenomenon  is  more  realistically 
represented,  with  particular  emphasis  on  higher  frequencies,  and  to  compare 
results  with  those  from  a  representation  based  on  St.  Venant's  theory. 
Similiarly,  Reference  8  examined  the  effects  of  the  constraint  of  waroing  on 
the  flutter  for  these  models  of  composite  wings. 
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Ouring  the  discussion  of  the  first  quarterly  Technical  Report  (7,8),  which  took 
place  on  5  September  1985  at  the  office  of  Dr.  Anthony  K.  Amos  (the  program  mon¬ 
itor  and  a  Program  Manager  for  AFOSR),  questions  that  were  raised  included  one 
that  dealt  with  the  following  observation:  in  studying  the  effects  of  the 
warping  constraint  on  the  dynamics  and  aeroelastic  stability  of  plate-like 
models  of  composite  aircraft  wings,  using  the  new  methodology  developed  at 
Fairchild  Republic  Company  (FRC),  along  with  wing  displacement  assumptions  in 
which  chordwise  curvatures  were  ignored,  it  appeared  startling  to  notice  from 
the  results  obtained  that  the  dynamic  and  aeroelastic  characteristics  of  these 
wings  are  dependent  on  their  chordwise  stiffnesses. 

During  the  month  of  September,  a  more  detailed  explanation  of  this  apparent  con¬ 
tradiction  was  formulated  and  the  following  conclusions  were  arrived  at: 

a)  The  dynamic  and  aeroelastic  characteristics  of  these  aircraft  wings  are 
only  implicitly  (not  explicitly)  dependent  on  the  chordwise  stiffnesses. 

b)  This  implicit  relationship  is  brought  about  by  Poisson  effects. 

c)  The  new  methodology  appears  to  expose  certain  details  that  may  not  be  vis¬ 
ible  when  conventional  methods  are  employed  to  analyze  composite  aircraft 
dynamic  and  aeroelastic  phenomena. 

The  difficult  task  of  determining  closed-form  solutions  to  the  dynamic  problem 
of  composite  aircraft  wings  with  the  coupling  and  warping  constraints  was  also 
begun  this  quarter.  Encouraging  progress  has  been  achieved  on  this  task,  which 
apparently  has  never  been  previously  accomplished.  The  accomplishment  of  the 
task  will  provide  the  aerospace  scientific  community  with  new  physical  insights 
into  these  dynamic  and  aeroelastic  phenomena.  It  should  also  provide  useful 
background  for  properly  investigating  the  effects  of  chordwise  curvatures  on 
these  phenomena,  which  is  planned  to  be  the  next  task. 
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5.4.2  Formulation 

Consider  an  aircraft  wing  fabricated  of  composite  materials  and  mathematically 
idealized  as  a  cantilevered  plate  subjected  to  an  aerodynamic  flow  over  its  sur¬ 
faces.  The  mathematical  statement  of  the  virtual  work  theorem  for  such  a  plate 
model  is  well  known  and  documented.  It  is  also  known  that  such  mathematical 
statements  of  the  virtual  work  theorem  for  a  laminated  plate  model  are  char¬ 
acterized  with  the  existence  of  so  many  variables  (in  the  statement), 
reflecting  the  various  directional  properties  for  the  laminated  plate  model, 
which  would  tend  to  interfere  with  any  physical  insight  that  might  be  desired 
from  a  phenomenological  analysis  employing  such  a  mathematical  statement.  The 
newly  discovered  affine  transformation  concept  (5,6)  was  developed  principally 
to  resolve  such  a  problem. 

This  new  concept  therefore  can  be  used  to  evolve  the  mathematical  statement  of 
the  virtual  work  theorem  in  an  affine  space  given  by  the  following  equation. 


iUQ  =  0  *  2 


L'SI  I'- 


X  X 
0  0 


)c  ♦  2D* 


(i  -  «)(w,x  v  Y 
Vo 


(l) 


2 

+  €w’x  X  w’y  y  1  +  (w’y  y  }  +  Llw’x  x  wx  y  +  L2*’v  y  w’x  v  t  <*x  dy  dt 

o  o  Voj  Vo  Vo  Vo  2  Vo  \yc  ]  o  yc 

■  if* jj  4  J J J 


OYIBO: 2:8/8/86 
VT3/ENG:17 


'U  '  “12 


^  D1 1°22 ' 


(DnD22)- 


Ll  *  f  n  ^3/4fp  ,1/4 
(D11}  (c22; 


(D11)1/4(D22)3/4 


•  O  * 
1  *0 


Dij  are  the  elastic  constants,  p  ,  is  the  material  density,  ^p  is  the  differ¬ 
ential  pressure  distribution,  w  is  the  displacement,  t  is  the  time,  A  integrals 
represent  area  integrals  and  h  is  the  wing  box  depth. 


Equations  (1)  and  (2)  therefore  form  the  basis  of  the  newly  developed  method¬ 
ology.  The  equations  of  motion  of  a  plate  model  of  an  aircraft  wing  can  now  be 
derived  by  prescribing  a  realistic  wing  displacement  and  using  Equation  (1). 


When  Equation  (1)  is  compared  to  its  physical  space  counterpart,  it  is  seen 
that  Equation  (1)  has  fewer  variables.  It  is  also  seen  that  Equation  (1)  con¬ 
tains  only  non-dimensionalized  stiffness  quantities  (compared  to  dimensional 
stiffness  quantities  in  its  physical  space  counterpart).  Another  feature  of 
this  new  methodology  which  makes  it  unique  is  that  the  non-dimensionalization 
(a  consequence  of  the  affine  transformation)  is  accomplished  before  assuming 
the  wing  deformations.  This  means  that  the  non-dimensionalization  is  independ¬ 
ent  of  how  the  wing  deforms.  A  non-dimensionalization  scheme  that  depends  on 
a  particular  assumption  of  the  wings  deformations  could  lead  the  analyst  to  an 
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incorrect  physical  interpretation  of  results,  since  the  wing's  deformations 
assumptions  have  inherent  errors  because  they  are  based  on  the  analyst's  judg¬ 
ments  and  experience.  This  observation  may  become  clearer  during  the  evolution 
of  a  warping  parameter  with  which  to  study  the  effects  of  the  warping  con¬ 
straint  phenomenon  on  the  status  and  dynamics  of  a  wing  fabricated  of  composite 
materials  later  in  this  study. 

If  the  chordwise  curvature  is  neglected  in  an  initial  approximation,  the  wing's 
deformation  may  be  assumed  as  follows: 

w(t'xo'yo)  *  ho(t’V  +  xoao{t‘yo)  (3^ 

where  h0  and  a  0  are  the  bending  and  twisting  displacements  respectively. 

It  can  be  shown  that  when  Equation  (3)  is  substituted  into  Equation  (1)  and  the 
variational  calculus  is  carried  out  for  arbitrary  hQ  and  &0,  the  following 
equations  of  -motion  are  obtained: 
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5.4.3  Evolution  of  Warping  Parameters 

The  evolution  of  the  warping  parameter  with  which  to  study  the  aeroelastic  warp 
ing  constraint  phenomenon  for  wings  fabricated  of  composite  materials  is  a  pro¬ 
cess  that  depends  on  the  sophistication  of  the  wing's  mathematical  model; 
whether  coupling  effects  are  included,  whether  the  wing's  chordwise  curvatures 
are  included  and  so  on.  Therefore,  any  warping  parameter  is  as  good  as  the  cor 
responding  wing's  displacements  assumptions.  However,  Equation  (1)  makes  it 
possible  for  the  analyst  to  determine  its  effective  independent  variables  even 
before  the  displacement  assumptions  are  made. 
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By  non-dimensionalizing  the  spanwise  space  variable  in  Equation  (4),  depending 
on  whether  one  is  interested  in  the  static,  dynamic,  coupled  or  uncoupled  dis¬ 
placements,  one  of  the  following  warping  parameters  may  be  useful. 
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q/ Cq)  is  defined  as  the  wing's  effective  aspect  ratio  and  D0*  and  L  are  the 
generalized  stiffness  and  coupling  ratios  respectively  (defined  in  earlier 
work  such  as  References  5  and  6) . 

Equations  (7)  thru  (9)  represent  the  appropriate  warping  parameter  for  dynamic 
deformation,  static  displacement  with  elastic  cross-coupling,  and  static  defor¬ 
mation  with  "geometric"  coupling  (  e  #  1)  respectively. 
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It  was  discovered  in  References  7  and  8  that  evolving  the  warping  parameter  in 
a  manner  shown  in  Equations  (1)  thru  (3),  should  enable  one  to  effectively  in¬ 
vestigate  the  effects  of  warping  on  the  composite  wing's  dynamics  (or  the  accu¬ 
racy  of  St.  Venant's  theory.  From  the  lamination  theory  for  composites  it  is 
known  that  while  D0*  and  (, lQ/zQ)  are  always  positive,  L  and  T2  can  be  positive 
or  negative.  However,  from  Equations  (1)  and  (3),  it  is  clear  that  whether  a 
composite  wing  has  positive  or  negative  coupling,  the  warping  effect  (in  terms 
of }  c)  is  unchanged. 

5.4.4  Results  and  Conclusions 

By  using  the  evolved  warping  parameters  defined  in  Equations  7  and  8  and  appro¬ 
priate  boundary  conditions,  the  boundary  value  or  the  eigenvalue  problems  asso¬ 
ciated  with  Equation  4  are  solved  in  a  closed-form  manner  to  determine  the 
wing's  coupled  vibration  frequencies.  The  presence  of  the  elastic  coupling  was 
found  to  significantly  complicate  the  mathematical  analysis,  required  to  solve 
the  problem  in  a  closed-form  manner.  This  fact  therefore,  is  probably  the  prin¬ 
cipal  reason  why  no  one  seems  to  have  been  successful  previously  in  obtaining 
closed-form  solutions.  The  necessity  for  closed-form  solutions  stems  from  what 
appears  to  be  a  lack  of  identified  coherent  trends  which  tend  to  characterize 
the  results  obtained  from  most  of  the  numerical  methods.  That  is,  the  earlier 
numerical  solutions  seem  to  be  a  collection  of  individual  isolated  results  with 
no  discernible  trends  exhibited  which  would  unify  the  meaning  of  the  results. 

In  the  process  of  trying  to  accomplish  this  task,  the  following  phenomena  were 
discovered:  (a)  the  coupling  mechanism  transforms  the  bending-torsion  oscilla¬ 
tions  (which  are  separately  conservative  problems)  into  a  non-conservative  prob¬ 
lem,  resulting  in  a  complex  determinant  for  the  eigenvalues;  (b)  one  "search" 
variable  (rather  than  two)  can  be  used  to  determine  the  eigenvalues  for  zeros) 
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of  this  complex  determinant,  and  (c)  the  frequency  parameters  are  only  func¬ 
tions  of  an  effective  aspect  ratio  (£0/c0),  a  nondimensional ized  coupling  param¬ 
eter,  12,  and  a  non-dimensional ized  twisting  stiffness  parameter,  Do*,  which  is 
proportional  to  the  warping  parameter  derived  earlier. 

Figure  1  shows  a  typical  behavior  of  the  real  and  imaginary  parts  of  the 
determinant,  which  coalesce  at  its  eigenvalues  (zeros). 

These  figures  therefore  depict  the  mechanism  by  which  coupling  (L2)  and  the 
warping  constraint  (proportional  to  0*o)  transform  the  first  vibrational  mode. 
The  word  transform  has  been  utilized  in  order  to  make  the  point  that  the  physi¬ 
cal  mechanism  known  as  coupling  transforms  the  nature  of  the  first  vibrational 
mode  in  the  process  of  transfering  energy  to  other  modes.  The  following  obser¬ 
vations  can  be  obtained  from  Figures  2-7:  the  first  coupled  frequency  parame¬ 
ter  (i)  increases  with  decreasing  coupling  (L2),  (ii)  increases  with  increasing 
Do*,  (iii)  is  independent  of  Do*  for  very  small  L2,  (iv)  increases  more  rapidly 
with  Do*  for  larger  coupling  (v)  decreases  more  rapidly  with  increasing  I2  for 
lower  Do*.  Observation  (i)  is  supported  by  data  computed  numerically  at  MIT4, 
observation  (ii)  implies  that  the  first  coupled  frequency  which  is  supposed  to 
be  dominated  by  the  first  bending  frequency,  becomes  more  "torsional"  as 
coupling  is  increased,  observation  (iii)  indicates  that  as  I2  approaches  zero, 
the  first  coupled  frequency  approaches  the  pure  (uncoupled)  bending  frequency, 
hence  it  becomes  independent  of  Do*;  this  might  provide  the  basis  to  explain 
the  modal  "transformations"  reported  for  instance,  by  the  study  at  MIT9. 

Figures  2-7  therefore  seem  to  suggest  a  potentially  accurate  method  to  trace 
the  origin  of  a  coupled  mode. 
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5.5  AEROELASTIC  DIVERGENCE  OF  ORTHOTROPIC  PLATES  by  Terrence  A.  Weisshaar 
5.5.1  Introduction 

The  aeroelastic  stability  of  lifting  surfaces  constructed  of  laminated  com¬ 
posite  materials  has  been  the  subject  of  numerous  research  studies  during  the 
past  decade.  A  variety  of  mathematical  models  have  been  used  for  these 
studies;  the  choice  of  model  has  been  dependent  upon  the  objectives  of  the 
researcher.  For  instance,  for  preliminary  design  or  trend  studies  the  level 
of  complexity  is  usually  kept  to  a  minimum,  consistent  with  the  experience  of 
the  analyst.  Aeroelasticians  have  only  recently  begun  to  realize  that  some  of 
their  lengthy  experience  or  "rules  of  thumb"  accumulated  with  metallic,  iso¬ 
tropic  structures  may  not  be  strictly  applicable  to  advanced  composite  mater¬ 
ials.  In  particular,  shearing  strains  in  laminated  materials,  while  safely 
neglected  for  certain  loadings  of  metallics,  assume  great  importance  in  com¬ 
posites.  For  instance,  the  ratio  E/G  is,  for  a  homogeneous  isotropic  mater¬ 
ial,  fixed  t'  be  2(1  +  H  )  =  2.6;  a  similar  ratio  for  a  laminated  material  may 
be  of  the  order  of  20  to  50.  As  such,  effects  of  shearing  are  likely  to  prop¬ 
agate  quite  far  from  their  source.  This  "source"  is  usually  a  restrained 
boundary  or  a  concentrated  load. 

To  Illustrate  the  potential  importance  of  shear  flexibility,  this  study  will 
examine  the  idealized  problem  of  static  aeroelastic  divergence  of  an  unswept, 
constant  property,  rectangular  flat  plate  constructed  of  laminated  composite 
material.  The  laminate  Is  assumed  to  be  orthotropic  with  respect  to  the  x-y 
axis  system  shown  In  Figure  1. 

The  stiffness  characteristics  of  the  plate  are  given  by  the  plate  stiffness 
matrix  [D^].  The  out-of -plane  deformation  w(x,y)  is  positive  downward. 
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along  the  +z  axis  formed  by  applying  the  right-hand  rule  to  the  x-y  axes.  A 
moving  airstream  with  velocity  V  indicated  in  Figure  1  creates  lift  that  acts 
on  this  plate. 

The  divergence  of  such  a  structure  has  been  studied  extensively.  The  most 
common  academic  problem  models  the  plate  as  a  beam-like  structure  using  Engi¬ 
neers'  beam  bending  and  twist  theories  (Bredt-Batho  theory).  This  theory 
neglects  the  inplane  (x-y)  normal  stresses  and  restraint  caused  by  the  fixed 
edge  of  the  plate  and  assumes  that  the  cross-sections  are  free  to  warp  out  of 
the  y-z  plane.  Neglect  of  this  so-called  “restraint  of  warping"  effect  does 
not  lead  to  serious  consequences  for  metallics  if  the  planform  aspect  ratio  is 
greater  than  about  6.  This  rule  of  thumb  (aspect  ratio  greater  than  6) 
appears  to  originate  from  an  early  study  by  Reissner  and  Stein  [1],  who  inves¬ 
tigated  the  effect  of  “restraint  of  warping"  upon  plates  with  a  variety  of 
loadings,  none  of  them  aeroelastic.  Whether  or  not  this  rule  applies  to  lami¬ 
nated  structures  is  of  concern  to  modern  structures  researchers. 

Crawley  and  Dugundji  [2]  as  well  as  Oyibo  and  Berman  (3]  shared  this  concern 
about  the  effect  of  material  properties  upon  the  restraint  of  warping  effect 
in  flat  plates.  They  studied  the  effect  of  a  warping  parameter  upon  the 
torsional  natural  frequencies  of  a  rectangular  flat  plate.  They  concluded 
that,  in  addition  to  planform  aspect  ratio,  the  plate  stiffness  ratio 
D,./D,,  must  be  considered.  This  ratio,  for  isotropic  metallic  materials, 

I  I  00 
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However,  for  a  laminated  plate,  It  is 


'11  ill 


Qn  y  1 


66 


.  S  .  q66  *1 

1*1 


A 


i  + 1 


z2dz 


(2) 


(3) 


The  integral  in  Eqn.  3  extends  over  the  thickness  of  the  ith  ply  comprising  an 
element  of  the  laminate.  The  quantities  and  are  ply  elastic  prop¬ 
erties  measured  along  orthotropic  directions.  These  properties  vary  with  ply 
orientation  with  respect  to  the  x-y  axis  system  in  Figure  1.  The  ratio  in 
Eqn.  2  may  be  of  far  different  magnitude  than  that  in  Eqn.  1. 


Reference  2  solved,  in  approximate  form,  the  equations  for  free  torsional 
vibration  of  a  laminated,  orthotropic,  cantilever  plate.  Included  in  the 
analysis  was  the  definition  of  a  non-dimensional  warping  parameter,  fi  . 
References  2  and  3  show  that  the  Inclusion  of  restraint  of  warping  effects 
leads  to  significantly  higher  torsional  frequencies  as  compared  to  those  pre¬ 
dicted  using  Bredt-Batho  or  St.  Venant  torsion  theory. 


5.5.2  Objective 

To  illustrate  the  potential  importance  of  laminate  geometry  to  even  the  sim¬ 
plest  aeroelastic  stability  problem,  an  example  was  studied.  Results  of  this 
study  were  then  compared  to  results  using  an  elementary  classical  solution. 
Potentially  Important  features  of  excluding  “restraint  of  warping"  effects  are 


1 1  lustrated. 
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The  example  involves  the  prediction  of  the  divergence  dynamic  pressure  of  an 
unswept,  rectangular,  orthotropic  laminated  plate.  An  idealization  using 
plate  theory  to  model  the  structural  deformation  was  used,  together  with  aero¬ 
dynamic  strip  theory  to  model  the  airloads  resulting  from  deformation.  A 
closed-form  solution  for  the  divergence  dynamic  pressure  is  found  and  then 
compared  to  that  obtained  using  conventional  St.  Venant  torsion  theory. 

5.5.3  Analysis 

The  variational  calculus  is  used  to  generate  the  governing  ordinary  differen¬ 
tial  equations  (in  terms  of  x)  and  associated  boundary  conditions.  The  fol¬ 
lowing  definitions  are  used: 


K  =  2c0 


16 


El  =•  cO 
GJ 


11 
4c0 
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°iic' 
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The  governing  differential  equations  are: 


EIw™  +  K01 


♦  qcaQ0  *  0 


So* 


IV 


GJ0 


II  III 

-  *wo 


qceaQ0  =  0 


The  associated  boundary  conditions  are: 

dw 


At 


0:  5^  -  0;  w0  -  0;  9  -  0;  $  - 


=  0 


(4) 

(5) 

(6) 


The  last  boundary  condition  at  x  =  0,  0'  (0)  =  0,  is  the  result  of  the  require- 
dw 

ment  that  _ o  (0)  =  0  at  all  points  along  the  clamped  edge.  This  produces  two 

dx 

BC's  for  a  plate  Instead  of  the  usual  single  beam  clamping  condition.  Boundary 


conditions  at  x  =  L  are: 
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EIw“  +  K©'  =  0 

111  II 

£Iwq  +  K0  =  0 

sV'=o  <7> 

o 

III  I  II 

-S0  +  GJ0  +  Kw„  =0 

0  0 

These  equations  are  greatly  simplified  if  the  bend/twist  stiffness  cross¬ 
coupling  parameter,  K,  is  zero,  implying  that  the  plate  is  orthotropic  with 
respect  to  the  x-y  axis  system.  With  K  =  0,  the  diverence  equations  for  0  can 
be  solved  independently  of  those  for  wQ.  The  $  equations  are: 


S  0IV  -  GJ#'  '  -  qcea  0  =  0 
o  o 


(8) 


0(0)  -  0 ' (0)  =  0 


(9) 


•  "(L)  =  So0,,'(L)  -  GJ0(L)  =  0  (10) 

The  parameter  SQ  is  indicative  of  the  retention  of  the  restraint  of  warping 

effect  in  this  one-dimensional  analysis.  If  SQ  is  zero  then  Eqn.  8  reduces 

to  the  classical  differential  equation  for  torsional  divergence  using  St. 

I 

Venant  torsion  theory.  In  addition,  if  SQ  is  zero,  then  the  BC's  Q  (0)  =• 

I  I 

0  and  0  (L)  =  0  are  not  present. 


To  solve  the  equation  set,  Eqns.  8  -  10,  we  first  non-dimensionalize  by  defin¬ 
ing  a  new  variable,  n,  as  follows: 


rj 


(ID 
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In  addition,  an  eigenvalue  X  is  defined  as 

X  2  -  qceaQL2/GJ  (12) 

A  nondimensional  warping  parameter  used  in  Reference  2  is  defined  as: 


(13) 


This  non-dimensional  parameter,  fi,  is  defined  in  terms  of  plate  constants  and 
geometry  as 


o-  A  /Ullwci 
P  48  'n  J'  J' 


u66 


(14) 


The  governing  differential  equation  is  now  written  as: 

-e"  -  2e  *  o  os) 

with  boundary  conditions: 

0  =  e'  =  o  at  rj  =  0  (16) 

and 

$"  =  09  "  -  9  =  0  (17) 

The  derivatives  in  these  latter  three  equations  are  with  respect  to  rj. 


Using  standard  techniques,  it  can  be  shown  that  the  characteristic  equation 
satisfying  Eqns  (15  -  17)  for  divergence  is 

2  +  (~  -  f)sin(g)sinh(f)  +  (  \  \  )cos(g)cosh(f )  =  0  (18) 

9  r  g  f 
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Where  f,  g  and  are  related  as  follows: 


:2  =  g2  +  i /fi 

(19) 

A2  =  g2d  +  fig2) 

(20) 

A2  =  f2(0f2  -  1) 

(21) 

Solution  for  A  proceeds  as  follows.  First,  choose  a  value  of  fi.  Because 
of  Eqn.  19,  with  fi  fixed,  only  g  or  f  is  now  independent.  Choose  g  as  the  in¬ 
dependent  parameter  and  iterate  to  solve  Eqn.  18.  This  was  accomplished  using 
a  short  computer  code.  A  plot  of  A  (divided  by  XQ)  versus  fi  is  provided  in 

i  5  2 

Figure  2.  The  parameter  is  equal  to  ir  /4  and  corresponds  to  the 
eigenvalue  obtained  when  fi  =  0,  the  classical  torsional  divergence  case 
with  St.  Venant  torsion  theory. 

2 

The  dynamic  pressure  at  divergence  is  proportional  to  X  •  Figure  2  indi¬ 
cates  that  the  inclusion  of  the  restraint  of  warping  effect  increases  the 
divergence  dynamic  pressure.  More  importantly,  failure  to  account  for  the 
restraint  of  warping  effect  most  likely  will  lead  to  predicted  divergence 
dynamic  pressures  that  are  lower  than  those  to  be  found  in  experimental 
tests.  While  Figure  2  indicates  the  effect  of  fi  upon  the  prediction  of  diver¬ 
gence  dynamic  pressure,  it  does  not  provide  adequate  insight  into  the  physical 
behavior  of  the  system.  To  acquire  a  more  complete  understanding  of  this 
restraint  of  warping  effect  on  plate  divergence,  let  us  consider  an  additional 
set  of  examples. 
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Crawley  and  Dugundji  [2]  also  present  experimental  data  related  to  free 
vibrations  of  flat  orthotropic  plates.  Their  specimens  include  8  plies  thick 
Hercules  ASI/3501-6  graphite  epoxy.  Plate  stiffness  (D.^)  data  from  Ref.  2 
can  be  used  to  illustrate  the  importance  of  orthot-opicity  to  restraint  of 
warping  effects. 

Table  1  shows  0^,  D&6  data  for  3  laminates.  The  first  laminate  is 

designed  to  be  stiff  in  bending,  the  second  laminate  is  quasi-isotropic  and 

the  third  laminate  is  stiff  in  torsion.  Plate  aspect  ratio,  L/c,  is  equal  to 

3.  This  data  was  used  to  compute  parameters;  eigenvalues  were  then 

calculated.  The  values  A. /X  are  shown  in  the  fourth  column  of  Table  1. 

o 

TABLE  1.  LAMINATE  STIFFNESS/DIVERGENCE  CHARACTERISTICS 


LAMINATE 

On /066 

£ 

2  ^/-P- 

[02/±30°]s 

16.94 

0.0392 

1.265 

[0/+45/90]s 

5.59 

0.0129 

1.14 

[i45/+45]s 

1  .216 

0.0028 

1  .059 

The  final  column  of  Table  1  shows  the  ratio  of  the  eigenvalue  X  (divergence 
speed  is  proportional  to  X  )  divided  by^/2.  This  value,  TT/2,  is  the  lowest 
value  of  X  that  would  be  computed  if  SQ  (or/S  )  were  zero.  Thus,  for  these 
three  laminates,  the  inclusion  of  plate-like  restraint  of  warping  effects  re¬ 
sults  in  higher  airspeeds  being  predicted  for  divergence  than  if  {}  were 
zero.  For  the  three  laminates,  the  two  predictions  (with  L/c  =  3)  differ  by 
26.5%,  14%  and  6.9%,  respectively.  To  bring  the  26.5%  error  down  to  6.9%,  the 
first  laminate  would  have  to  have  its  aspect  ratio  increased  from  3  to  10.8! 
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The  mode  shape  at  divergence  Is 


8Cn)  =  A] (sin(gTj)  -  fsinh(ft))) 
+  A2(cos(g^)  -  cosh( fTj) ) 


(22) 


From  boundary  conditions,  the  constant  Aj  is: 


-A1(g2sin(g)  +  gfsinh(f)) 

A2  - - 

2  2 
(g  cosg  +  f  coshf) 


i(n) 


e<i\) 

eo7 


(23) 

(24) 


If  restraint  of  warping  effects  are  excluded,  this  mode  shape  would  simply  be 


0(T))  =  Asin(^T^)  (25) 

The  mode  shapes  or  torsional  deformation  at  divergence  are  both  interesting 

and  informative.  Figure  3  shows  the  normalized  mode  shape  given  in  Eqns.  22  - 

24  for  several  values  of  ft  .  Large  values  of  ft  are  associated  either  with 

low  aspect  ratio  plates  or  with  large  values  of  D  /D  .  Also  shown  is 

11  66 

the  classical  St.  Venant  theory  divergence  mode  shape  (  ft-  0).  Even  for 
small  values  of  ft  ,  the  mode  shape  at  divergence  is  significantly  different 
than  the  quarter-sine  wave  associated  with  St.  Venant  torsion  theory.  Differ¬ 
ences  are  particularly  noticeable  In  the  region  near  t)  *  0.  The  idealization 
used  In  this  study  predicts  airloads  that  are  proportional  to  the  mode  shape 
for  0(T))  plotted  in  Figure  3.  The  Increased  divergence  dynamic  pressure  that 
results  when  restraint  of  warping  is  considered  Is  seen  to  be  due  to  reduced 
Incidence  In  the  Inboard  region.  When  equals  unity  there  is  a  reversal  of 
incidence  in  the  region  0  <  n  <  0.75  which  leads  to  very  high  divergence 
airspeeds. 
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5.5.4  Conclusions 

This  study  has  examined  the  aeroelastic  divergence  of  flat  cantilever  ortho¬ 
tropic  plates.  These  results  are  compared  to  those  obtained  using  classical 
St.  Venant  torsion  theory  that  neglects  plate  bending  terms  that  reflect  re¬ 
straint  of  warping  of  plate  cross-sections.  A  non-dimensional  parameter  £  is 
defined  and  used  to  determine  the  degree  to  which  warping  effects  are  impor¬ 
tant.  This  parameter  is  composed  of  the  ratio  of  spanwise  bending  stiffness 
to  torsional  stiffness  and  plate  structural  aspect  ratio.  For  values  of  £> 
0.01,  the  inclusion  of  warping  restraint  is  found  to  result  in  larger  dynamic 
pressures  at  divergence  than  would  be  found  using  St.  Venant  torsion  theory. 
Thus  St.  Venant  torsion  theory  leads  to  overly  conservative  results.  Equally 
important,  for  £>0.01  the  divergence  mode  shape  with  restraint  of  warping 
included  may  be  substantially  different  from  that  found  using  St.  Venant  tor¬ 
sion  theory.  This  latter  observation  is  important  to  the  selection  of  assumed 
modes  for  approximate  solutions  to  aeroelastic  stability  problems. 


The  differences  encountered  In  this  divergence  study  are  due  to  the  ortho¬ 
tropic  behavior  of  plates  and  are  not  confined  to  laminated  composites.  How¬ 
ever,  some  laminate  designs  display  an  additional  feature  known  as  bend/twist 
stiffness  cross-coupling.  This  effect  couples  together  the  twisting  and  out- 
of -plane  bending  of  the  plate  in  such  a  fashion  that  the  full  set  of  equa¬ 
tions,  described  earlier,  must  be  solved.  A  second  limitation  of  the  present 
results  is  that  camber  bending  has  been  suppressed  by  truncation  of  the  series 
given  in  Eqn.  8.  Inclusion  of  additional  terms  may  be  necessary  to  further 
refine  the  solution. 
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5.6  A  STUDY  OF  THE  EFFECTS  OF  PLATE  STRUCTURAL  IDEALIZATIONS  UPON 
TAILORING  METHODOLOGY  by  Terrence  A.  Weisshaar 

The  great  wealth  of  analytical  studies  related  to  aeroelastic 
tailoring  that  have  appeared  recently  have  added  insight  and  confidence 
to  the  advanced  composite  design  integration  process.  Many  different 
structural  idealizations  are  used  to  demonstrate  tailoring  features  and 
potential  improvements  to  be  found  through  tailoring.  Some  of  these 
idealizations,  while  accurate  for  conventional  metallics,  have  been 
found  to  be  lacking  in  essential  features  when  advanced  composite  con¬ 
struction  is  the  mode  of  construction.  Elimination  or  disregard  of  some 
elastic  effects  that  may  safely  be  neglected  for  metallics  has  been 
found  to  be  unwise  in  the  study  of  tailored  structures. 

The  study  of  the  effects  of  tailoring  as  it  affects  structural  or 
aeroelastic  dynamic  response  is  conducted  along  two  parallel  approaches. 
The  first  approach  is  to  select  a  certain  geometric  design,  for  instance 
a  fixed  planform,  and  to  choose  a  specific  material,  e.g. 

G raphite/ epoxy .  The  tailoring  study  then  examines  the  effects  of  ply 
orientation,  thickness  and  stacking  sequence  upon  the  response  of 
interest.  The  results  are  valid  only  for  that  material  and  geometry  and 
stacking  sequence/ thickness  combination. 

A  second  approach  is  to  define  a  set  of  "natural  variables"  that 
are  independent  of  any  actual  planform,  material  or  geometry.  Studies 
then  identify  desirable  combinations  of  these  natural  variables  and 
ass' ss  their  effects  upon  response.  Such  an  approach  has  advantages  and 
disadvantages.  A  major  advantage  is  that  special  behavior  can  be  attri¬ 
buted  to  a  distinct  generic  feature  of  laminated  composites.  When  lam- 
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inate  parameters  are  used,  for  instance,  in  a  dimensional  analysis  that 
uses  fiber  sweep  as  a  variable,  a  number  of  important  stiffness  vari¬ 
ables  change  simultaneously.  One  is  hard  pressed  to  attribute  desirable 
(or  undesirable)  changes  to  a  specific  cause.  The  use  of  natural  vari¬ 
ables  that  are  nondimens ional  in  nature  eliminates  this  difficulty. 

There  are  two,  in  particular,  problems  with  the  use  of  natural 
variables.  First  of  all,  how  many  truly  independent  nondimens ional 
parameters  are  present?  Secondly,  having  identified  these  parameters, 
what  are  their  true  physical  limits  or  ranges  of  value?  Both  questions 
are  important  since  any  analytical  model  consists  only  of  numbers  and, 
if  used  incorrectly,  may  model  an  impossible  physical  process.  A  case 
in  point  is  the  insertion  of  a  negative  stiffness  into  a  vibration  equa¬ 
tion  or  use  of  a  negative  viscous  damping  constant. 

This  study  exams  the  problems  associated  with  nondimens ional  model¬ 
ing  of  plate  vibration  idealizations.  The  first  section  examines  the 
circumstances  under  which  the  elements  of  the  matrix  of  plate 
stiffnesses  may  be  considered  to  be  independent.  The  second  section 
suggests  nondimens ional  parameters  for  use  in  these  types  of  scudies  and 
defines  their  number  and  limits.  The  third  section  illustrates  the  use 
of  these  parameters  in  a  free  vibration  study. 

5.6.1  Plate  Stiffness  Parameter  Independence  for  Tailorinq  Studies 

Lifting  surface  structural  idealizations  for  tailoring  studies  fre¬ 
quently  involve  the  use  of  equivalent  plate  assumptions.  In  this  case, 
there  are  six  plate  constants  (Djj,  D^,  D^,  D^,  D^,  D^)  to  be 
determined.  The  effects  of  changes  in  these  six  parameters  upon  free 


vibration  and  response  to  external  deformation  dependent  loads  is  of 
importance  to  design  optimization.  However,  before  any  parameter  varia¬ 
tion  can  be  studied,  it  must  be  determined  which  plate  stiffness  parame¬ 
ters  are  truly  independent  and  under  what  conditions  or  restrictions 
this  parametric  independence  arises.  An  examination  of  this  determina¬ 
tion  is  the  subject  of  this  section. 

To  become  familiar  with  the  most  basic  situation,  first  consider 
the  isotropic  plate  stiffnesses.  These  are: 


11 


16 


12 
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D22  -  D  -  Et3/12( 1  -  u2) 

D26-° 

UD 

D( 1  -  p)/2 


(1.1) 

(1.2) 

(1.3) 

(1.4) 


Isotropic  plate  design  uses  three  independent  design  variables,  E,  t  and 
u.  Material  constitutive  relationships  do  not  allow  and  D22  to  be 
changed  independently.  It  would  appear  at  first  glance  that  three  of 

the  plate  stiffnesses  could  be  treated  as  independent,  but  this  is  not 

3  - 

the  case  since  E  and  t  appear  as  a  unit.  Let  us  specify  that  ■  D 

and  ■  G,  where  D  and  G  are  fixed  values.  Solving  for  the  parameters 

E,  tu,  to  construct  such  a  plate,  we  find: 


48D66(1 


48G(l  -  -) 
D 


Dj2  Is  now  fixed  to  be: 


(1.5) 

(1.6) 


0 


12 


-  D(  1 


(1.7) 


Note  also  that  the  equation  for  u  has  limits  -1  <  u  <  0.5  for  isotropic 
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G 

materials  so  that  the  physical  limits  on  choices  of  —  are: 

D 

T  <  £  <  1  (1.8) 

D 

This  means  that  physically  we  cannot  construct  a  plate  with  stiffness 

ratio  D, ,/D  outside  these  limits  if  we  are  restricted  to  isotropic 
00  i  1 

materials.  A  question  remains,  if  there  are  three  design  parameters 
(E,t,u)  why  can't  we  arbitrarily  fix  three  plate  constants?  This  ques¬ 
tion  will  be  addressed  later  in  the  discussion. 

This  example  indicates  that  two  plate  stiffnesses  are  independent 
while  the  third  stiffness  depends  upon  the  other  two  because  of  material 
constitutive  relations.  We  also  see  that  there  is  only  one  independent 
ratio  of  stiffnesses,  that  is,  either 


22 

>1! 


1  -  2- 


_66 

’ll 


D66  °12 

or  - —  -  1  -  2- — 

D11  D11 

Next,  let  us  consider  an  orthotropic  material  for  which 


(1.9) 

(1.10) 


y  11 

(1.11) 

11 

12 

Q12c3 

(1.12) 

12 

12 

Q22C 

(1.13) 

22 

12 

V3 

(1.14) 

*66 

12 

16 

m 

°26  "° 

(1.15) 

Note  that  t  appears  in  every  parameter  and  it  cannot  be  used  as  an 
independent  design  parameter  to  adjust  relative  values  of  the  D^'s. 


1 

i 
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For  the  orthotropic  plate  there  are  four  material  parameters,  Q^,  Qj2» 
^22*  ^66’  chese  material  parameters  are  independent,  then  all 

four  plate  stiffnesses  may  be  prescribed  independently.  Consider  the 
case  for  which 


(1.16) 

(1.17) 

(1.18) 

(1.19) 


Let  Etl  *  E^t3,  E22  ■  E22t3‘  G12  *  **12^  wlth  ^  fixed*  If  Dn’  °22’ 

D._  and  are  the  prescribed  values  of  plate  stiffness,  then  the 
X  i,  00 

values  of  ^22,  u12  (or  u^)  and  necessary  to  guarantee  this 

selection  are: 


independent.  Notice  that  the  nondimens ional  ratio  Di2^11^22  aPPears  in 
two  of  these  equations.  If  we  were  to  require  values  of  Dj2,  D,  ^ ,  D^2 
that  led  to  values  of  this  ratio  in  excess  of  unity,  then  E^j  and  E22 
would  be  negative.  Such  cases  are  prohibited  by  thermodynamical  con¬ 
sideration,  namely  energy  conservation.  Additional  study  of  the  way  in 


f, 

/ 


which  the  plate  stiffness  parameters  may  be  varied  is  in  order. 

The  availability  of  laminated  composite  materials  introduces  design 
freedom  to  the  materials  selection  process.  Rather  chan  constructing  a 
material  with  engineering  constants  iC  is  more  *ikely 

that  a  Hooke's  Law  or  stiffness  matrix  [ Q ]  will  be  constructed.  Dis¬ 
cussion  of  this  construction  will  be  restricted  to  laminates  with  plies 
of  equal  thickness  and  a  fixed  overall  thickness. 

Let  us  first  consider  a  single  ply  of  thickness  t.  Place 
stiffnesses  are  defined  as: 


11  12 

,  .5a£ 
12  12 
-  3 


22  12 

,  .hd 

16  12 

,  .hd 

26  12 

,  .hd 

66  12 
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(1.25) 


(1.26) 


(1.27) 


(1.28) 


(1.29) 


The  bar  over  the  terms  indicates  that  they  have  been  transformed  by 
a  ply  angle  rotation.  The  quantities  are  defined  in  terms  of 
material  properties  and  angular  rotation  9  (Ref.  1). 


Qll  "  U1  +  U2  cos(20)  +  U3  cos(49> 
Q22  -  U  -  U2  cos(29)  +  U3  cos(40) 

QI2  -  U4  -  U3  cos(4  9) 

Q66  -  U$  -  U3  cos(4  9) 

Ql6-~i  U2  sin(2  0)  -  U3  sin(49) 


(1.30) 

(1.31) 

(1.32) 

(1.33) 


(1.34) 


(1.35) 


Q26  *  "  1^2  sin(29)  +  u3  sin(4  0) 

The  U1  terms  are  combinations  of  ply  material  properties, 
written  as: 

1 


They  are 


"l  'a0'3!!  *  *>22  *  *>12  *  “>66 > 
1 


U2  '  2«>ll  '  *>22  > 
1 


°3  '  8<Q11  +  »22  -  2<>12  -  “>66> 
0  -  * 


4  '  I^U  *  Q22  *  6<>12  '  4<W 
°S  •?«„  *  <=22  -  *>12  +  ‘V  -ra'l  -  V 


(1.36) 

(1.37) 

(1.38) 

(1.39) 

(1.40) 

Four  of  the  five  terms  are  Independent  and  invariant  for  a  specified 
material.  Material  selection  to  define  the  L^'s  is  analogous  to  choos¬ 
ing  E  and  u  for  an  isotropic  material.  This  material  choice  constitutes 
a  design  freedom.  Another  design  freedom  is  ply  orientation,  9.  Once 
the  single  material  is  prescribed,  for  a  laminated  material  only  four  of 
the  six  D^'s  may  be  chosen  independently.  Combining  equations  (1.24- 
1.40)  we  have: 


°ll  +  °22  +  20 12  -  2(U1  +  Vfl 


and  °66  -  °12 


<U5  -  Vfl 


(1.41) 

(1.42) 

Leaving  the  case  of  the  single  ply,  next  consider  a  laminate  with  4 
plies  of  a  given  material.  Restrict  the  discussion  to  symmetrical  lam¬ 
inates.  Tsai  and  Hahn  (Ref.  1)  show  that  the  plate  stiffnesses  may  be 
written  as: 


°11  "  VoD  +  VlD  +  V 


3D 


'12 


22 


U  V  -  U  V 
4V0D  U3V3D 


U1V0D  "  U2VID  +  U3V3D 


°16  *  -  2V2D  +  V4D 


(1.43) 

(1.44) 

(1.45) 

(1.46) 
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(1.47) 


D  «  -  -{J  V  +UV 
26  2^2  2D  3  4D 


66 


U5V0D  '  U3V3D 


(1.48) 


where  V, 


OD 

The  terms  V 


C3/12. 


ID’  V2D’  ^3D’  V4D  are  ^e^net*'  ^or  3  laminate  with  plies,  2n 


symme  trical , 


V  -  21  .ln<2B  )  s 
k*l 


V2D  '  22  Sl"<2\>  8k 

k“l 


V  -  21  cos(46  )  B 

k**l 


V  -  2  I  4l»(49k)St 
k-1 


(1.49) 

(1.50) 

(1.51) 

(1.52) 
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The  terms  0^  are  equal  to  -j{  t^  -  t^_j)  and  correspond  to  area  moments  of 


inertia  of  the  kth  plies  (1  <  k  <  n)  lying  between  positions  t^  ^  and  t^ 
from  the  plate  midsurface. 

3 


With  4  plies  (2  independent  layers),  a  given  material  and  with  t' 


fixed  (V0D  is  then  prescribed)  only  two  of  the  values  of  are 


independent.  This  is  so  because  Eqns.  1.41  and  1.42  are  valid  no  matter 

how  many  plies  are  used.  Let's  choose  D.,  and  D,,,  to  be  zero.  In  this 

lb  zb 


case , 


+  v  *  -U  V 
2T 2  2D  U3  4D 


(1.53) 


and 


-Lt  V 
2^2  2D 


U,  V . 

3  4D 


(1.54) 


The  only  solution  in  this  case  is  V ^  ■  V^D  -  0;  For  this  situation, 


0^  *  0,  i  90°,  corresponding  to  a  unidirectional  laminate  [0°,  0°]g  or  a 


cross-ply  laminate  [0°,  ±  90°]  .  The  parameters  D.  ,  ,  D, . ,  D.„,  and  D 

s  1  l  1 1  ll 
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are  now'  determined  since  no  other  design  variables  remain 


Adding  a  third  ply  (6  layers)  implies  that  three  values  are 

independently  "tailorable"  since  three  of  the  V  summations  in  Eqns. 

1.49-1.52  can  be  controlled.  For  instance,  we  might  choose 

*  0  and  ■  D^.  In  principle  three  values  of  6^  could  be 

found  to  create  this  situation  (note  that  sine  and  cosine  are  periodic). 

Extending  this  discussion,  four  plies  allow  the  choice  of  all  four 

values  of  V  and  the  four  independent  D  ,  can  be  specified. 

ij  s 

What  happens  when  the  number  of  plies  exceeds  four?  In  this  case, 
the  number  of  design  variables  exceeds  the  number  of  V  parameters. 
There  are  more  unknowns  (0^'s)  than  equations.  The  system  is  overdeter¬ 
mined  and  thus  the  0^  solution  set  to  fix  V^,  V^D ,  V^,  is  not 
unique.  This  is  one  great  advantage  of  advanced  composites,  since  there 
are  a  number  of  different  laminate  designs  with  equivalent  stiffnesses. 


5.6.2  Conclusion 

It  appears  that  the  six  plate  stiffness  parameters  (D^)  can 
regarded  as  independent  under  certain  conditions.  For  instance,  if  a 
certain  ply  material  is  chosen  then  four  or  more  plies  provide  tailoring 
independence  for  four  plate  stiffnesses,  within  the  limits  imposed  by 
material  relations.  Choosing  the  material  properties  then  adds  addi¬ 
tional  freedom  necessary  to  make  the  remaining  two  plate  stiffnesses 
independent.  If  the  parameters  are  to  be  regarded  as  independent 
then  one  must  allow  changes  in  material. 
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5.6.3  Appendix 


The  suppression  of  camber  bending  in  plate  vibrations  can  lead  to 
serious  errors  in  both  the  magnitude  of  the  natural  frequency  and  the 
mode  shape  determined  from  the  equilibrium  equation.  In  unpublished 
work  by  Liu  (not  related  to  Che  present  contract)  he  compares  the  plate 
frequency  with  the  results  using  the  classical  beam  bending  frequency 
expression  and  the  classical  beam  torsion  equation.  His  comparisons 
extend  over  a  wide  range  of  aspect  ratios,  but  are  restricted  to  ortho¬ 
tropic  plates.  A  summary  of  these  results  follows. 

A  rectangular,  cantilever  plate  with  the  following  dimensions  is 
considered.  The  chord  dimension  is  c*20"  while  the  span  dimension  is  L 
■  rc  *  20r  ,  where  r  is  the  plate  aspect  ratio.  Plate  stiffnesses  are 
fixed  to  be: 

°n  '  122-92  1Wn*  °16  *  °26  *  ° 

D12  -  33.66  lb-in.  D  -  46.5  lb-in. 

D22  -  801.59  lb-in. 

2 

The  mass  per  unit  area  is  p  ■  0.00171  slug/in  .  For  comparison,  bean 
parameters  are  defined  as: 


EI  ■  c(D22  -  D^/Dn) 
M  »  pc 


The  frequencies  of  the  first  two  bending  modes  are  predicted  to  be: 


flb  -  1319. 6/L2  and  f2  -  8269. 5/L2  (Hz.) 


For  torsional  vibration,  with  t  ■  plate  thickness, 


GJ  *  4cD66 

t  -  1  3/t  A  .  .2.  13 

J  T2PC  ^  +  )  *  f20c 

where  the  ratio  t/c  has  been  neglected.  The  first  two  torsional  fre¬ 
quencies  are  predicted  to  be: 

fu  -  49.47/L  (Hz.) 

f2t  "  148*4/L  (Hz‘) 

A  Rayleigh-Ritz  analysis  was  completed  for  12  different  values  of 
aspect  ratio,  r.  A  polynomial  series  representation  was  used  with  8 
terras  in  the  spanwise  direction.  In  the  chordwise  direction  the  number 
of  polynomial  terms  (denoted  as  n^)  was  varied.  If  n  •  1  then  chord- 
wise  rigidity  is  assumed.  If  n^  *  2  parabolic  camber  bending  is 
allowed.  Results  are  shown  on  the  following  pages.  In  each  case, 
columns  two  and  three  represent  the  ratio  of  natural  frequency  predicted 
from  Rayleigh-Ritz  plate  analysis  to  natural  frequency  predicted  from 
elementary  beam  analysis. 

These  results  illustrate  two  effects.  As  shown  by  previous 
authors,  torsional  frequency  error  is  strongly  dependent  upon  aspect 
ratio.  This  is  true  whether  camber  bending  is  neglected  or  not.  This 
difference  between  plate  results  and  beam  results  is  due  to  the  neglect 
of  inplane  warping  in  the  beam  analysis.  In  the  case  of  the  bending 
frequency,  as  aspect  ratio  increases,  the  discrepancy  between  plate 
theory  and  bending  theory  actually  increases  if  camber  bending  is 
suppressed.  The  discrepancy  is  diminished  with  increasing  aspect  ratio 
of  camber  bending  (n  ■  2)  is  permitted. 
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5.7  AN  EXTENDED  PLATE/BEAM  FINITE  ELEMENT  FOR  COMPOSITE  TAILORING  STUDIES 
by  Terrence  A.  Weisshaar 

There  is  a  growing  appreciation  of  the  importance  of  including 
inplane  shear  stiffness  effects  into  analytical  models  that  represent 
slender,  advanced  composite,  beam-like  structural  elements.  This  is  the 
so-called  "restraint  of  warping  effect"  discussed  in  a  number  of  recent 
studies  (1-3).  This  effect  stems  from  constraining  the  root  cross-section 
of  a  plate-like  or  beam-like  structure  to  remain  undeformed  in  its  own 
plane,  an  effect  that  is  usually  unaccounted  for  in  elementary  studies. 
Accountability  for  the  warping  constraint  may  be  handled  with  several 
levels  of  mathematical  complexity.  In  the  simplest  case,  one  might 
simply  resort  to  a  detailed  finite  element  model  of  an  anisotropic 
plate-like  lifting  surface.  Such  an  approach  is  cumbersome  and  may  be 
found,  after  great  effort,  to  be  unnecessary.  On  the  other  hand, 
restraint  of  warping  effects  are  automatically  Included  in  any  plate 
analysis . 

An  alternative  way  of  retaining  accuracy  while  maintaining 
computational  simplicity  is  to  develop  an  "extended  beam"  finite 
element.  By  extended,  it  is  meant  that  an  important  plate-like  feature, 
that  of  chordwise-camber,  is  retained  as  a  nodal  degree-of-freedom  to  be 
either  used  or  suppressed.  Such  elements,  with  limitations,  are 
discussed  by  Weisshaar  [4]  and  Chen  and  Yang  [5].  Weisshaar's  model  is 
deficient  in  that  it  does  not  include  transverse  shear  deformation,  nor 
does  it  include  sophislticated  shape  functions  such  as  those  used  by 
Chen  and  Yang.  Chen  and  Yang,  on  the  other  hand,  do  not  retain  camber 
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bending  as  a  degree  of  freedom.  As  a  result,  they  exclude  an  important 
effect  from  an  otherwise  useful  element. 

The  term  "extended  beam-element",  as  applied  to  composite  tailoring 
studies,  denotes  a  modified  beam  idealization  for  which  at  least 
parabolic,  chordise  camber  (bending)  is  retained.  Otherwise,  torsional 
stiffness  as  well  as  mode  shape  information  is  flawed.  As  a  result,  the 
accuracy  of  aeroelastic  stability  studies  may  be  compromised. 

This  appendix  outlines  the  basic  steps  necessary  to  develop  an 
extended  beam  finite  element  for  dynamical  studies.  Such  an  element  is 
essential  for  checking  the  results  of  other,  more  sophisticated,  studies 
as  well  as  generating  useful  trend  information. 

Consider  the  beam  element  shown  in  Figure  A-l.  The  element  has 
three  degrees  of  freedom  at  each  end.  These  are:  (1)  bending 
displacement;  (2)  bending  slope;  (3)  torsional  twist.  A  beam  cantilever 
condition  will  require  that  all  three  degrees  of  freedom  be  zero.  A 
plate  cantilever  condition  requires  that  the  edge  (perpendicular  to  the 
plane  of  the  paper)  remain  straight.  This  translates  into  the 
requirement  that  the  twist  rate  (first  derivative  of  torsional  twist)  be 
zero  and  that  no  camber  due  to  Poisson's  ratio  effects  be  present.  For 
the  model  shown  in  Figure  A-l,  such  constraints  are  uncontrollable. 

To  control  boundary  conditions  related  to  camber  and  warping 
restraints,  two  additional  degrees  of  freedom  at  each  node  are  required. 
One  degree  of  freedom  is  twist  rate  while  the  second  degree  of  freedom 
is  camber  curvature.  Let  us  turn  to  development  of  the  theoretical 
expressions  necessary  to  define  an  extended  beam  finite  element 


stiffness  matrix 


Shovm  in  Figure  A-2  is  an  element  of  length  L  with  five  degrees  of 


freedom  at  each  end.  These  are,  at  the  ends:  (1)  and  (6),  bending 
displacement;  (2)  and  (7),  bending  slope;  (3)  and  (8),  twist  rotation; 
(4)  and  (9),  twist  rate;  and,  (5)  and  (10),  camber  curvature.  An 
assumed  displacement  function,  w(x,y),  is  defined  to  be  (positive 
upward) 


w(x 


>y>  -  j*i(y)j  -  x|5i(y)j{qi}  +  x2{8i(y>}^i 


} 


(A.  1 ) 


In  Eqn.  A.l  the  functions  Y^,  6^,  8^  are  functions  of* the  spanwise 
coordinate,  y.  The  functions  Y^  provide  spanwise  displacement  due  only 
to  bending,  while  the  6  functions  give  information  on  twisting  effects. 
The  8^  functions  account  for  parabolic,  chordwise  camber  effects. 

Let  us  first  consider  a  beam  model  based  upon  laminated  plate 
theory.  In  the  case  of  a  mid-surface  symmetric  laminated  plate,  the 
relationship  between  plate  bending  moments,  twisting  moment  and 
curvatures  may  be  expressed  as  (see  for  instance,  Tsai  and  Hahn): 
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(A. 2) 


The  elements  are  anisotropic  plate  stiffness  constants  that  are 

functions  of  laminate  ply  geometry,  material  properties  and  stacking 

sequence,  while  <  ,  k  and  <  are  plate  curvatures, 
x’  y  xy  r 

If  we  adopt  the  y-axis  in  Figure  A-3  as  a  reference  axis,  the  plate 
deflection  w(x,y)  and  the  beam  deflection  h(y)  and  twist  a(y)  are 


defined  as: 


h(y)  -  w(0,y) 


(A. 3) 


and 


W  I  x-0  '  "a(y) 

In  this  case,  the  plate  curvatures  can  be  approximated  as  follows, 


(A. 4) 


and 


2  2 
<  *  3  w/3y 
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2  2 
■c  -  3Zw/3xZ 
x 

<  -2 
xy  3x3y 


The  strain  energy  density  of  the  plate  is  written  as 


Now,  let 


where 


K  <  < 

x  y  xy 


xy 


xy 


[A] {q} 


(A. 5) 

(A. 6) 
(A. 7) 


(A. 8) 


(A. 9) 


[A] 


28. 


/  #  » •  2  * ' 
Y.  “  x6.  +  x  6. 


+  2x6i 


and  the  notation  ()  denotes  differentiation  with  respect  to  y. 

As  a  result,  the  system  stiffness  matrix  is  given  by  the 
expression: 


(A. 10) 
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(A. 11) 


L  c/2  T 

I^i j]  -  /  /  [A]r[D^][A]dxdy 


o-c/2 


The  elemental  stiffness  matrix  depends  upon  the  choice  of  assumed 
displacement  shape  functions  Y^,  6^,  0^.  For  Y^>  choose 

Ti  -  1  -  3(f)2  -  2(f)3 
Y2  -  -y  +  2y2/L  -  y3/L2 

Y3  “  Y4  “  Y5  “  ° 

Y6  “  3(^)2  '  2(L)2 
Y?  -  -y3/L2  +  y2/L 

Y8  "  Y9  "  Y10  ”  0 

For  an  approximation  to  elemental  twist,  choose 


41  ■  52  *  0 
43'T1 

44  '  t2 

45  ■  0  ‘  S6 
48  '  \ 

4  9  •  h 


410  -  0 


For  camber  deformation,  let 


*1-02 


3  4 

S5  -  1  -  4(2)  +  3(2) 


S6  -  ®7 


S8  ■  S9 


8io  '  3(f>  -  2<f> 

After  multiplying  the  triple  matrix  product  shown  and  integrating 


with  respect  to  x,  we  have  the  following  relationship: 


The  elements  k 


fkijdy 

o  J 


ij* 


which  require  integration,  are  (note  k 


ij 
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An  examination  of  these  expressions  clearly  reveals  warping  and  other 

3  5 

edge  effects.  These  terms  are  multiplied  by  the  factor  c  / 1 2  or  c  /80. 
They  are  not  included  in  conventional  elements.  Note  that  the  term 
adds  torsional  stiffness  to  the  (3,3)  and  (8,8)  terms,  but  only  if  a 
nonlinear  torsional  deformation  function  is  used.  The  term  4cD,,  is  an 

DO 

equivalent  beam  GJ  stiffness  parameter. 
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6.0  OVERALL  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 

In  this  research  program,  the  phenomenon  of  constrained  warping  of  composite 
aircraft  wings,  resulting  from  aeroelastically-induced  loads,  has  been 
investigated.  In  the  preliminary  phase  (first  year)  of  this  research,  the  in¬ 
teraction  between  the  constrained  warping  and  the  statics,  free  vibration,  di¬ 
vergence  or  flutter  of  plate  models  of  composite  wings  without  elastic  coupling 
were  studied.  In  addition,  a  study  of  the  interaction  of  the  constrained 
warping  and  the  free  vibration  for  composite  wings  with  elastic  coupling  was 
initiated.  The  principal  results  revealed  by  these  preliminary  studies  include 
the  following,  (i)  the  constrained  warping  phenomenon  has  a  very  significant 
influence  on  the  static,  free  vibration  and  aeroelastic  stability  of  aircraft 
wings  fabricated  of  composite  materials,  (ii)  this  influence  may  be 
conservative  or  nonconservative  depending  on  the  configuration  of  the  wing, 

(iii)  there  exist  closed-form  solutions  to  the  problem  of  free  aeroelastic  os¬ 
cillations  of  composite  aircraft  wings  with  constrained  warping  and  elastic 
coupling,  (iv)  these  newly  discovered  closed-form  solutions  were  instrumental 
in  formulating  a  physical  explanation  (which  has  been  non-existent  until  now) 
of  the  vibrational  "modal  transformation"  phenomenon  in  models  of  composite  air¬ 
craft  wings. 

These  results  therefore  have  been  able  to  provide  answers  to  some  basic  ques¬ 
tions  about  the  fundamental  nature  of  the  warping  phenomenon.  However,  these 
new  results  have  raised  new  questions  as  well:  For  example,  does  the  existence 
of  closed  form  solutions  just  established  for  the  wing's  free  oscillations 
imply  the  existence  of  closed-form  solutions  for  the  divergence  or  flutter 
problems?  Other  questions  include:  is  it  possible  to  formulate  a  numerical 
super-element  which  can  accurately  capture  the  essence  of  these  dynamic 
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phenomena  which  can  be  used  for  analyzing  more  complex  models/configurations  of 

*  i 

composite  aircraft  wings?  j 

Some  of  the  issues  have  been  addressed  in  a  preliminary  fashion  during  the  cur¬ 
rent  phase  of  the  research  program.  In  a  future  research  program,  therefore, 
it  should  be  desirable  to  deal  with  these  issues  in  greater  depth  by,  (a) 
completing  the  analysis  of  the  elastically  coupled  free  vibration  with 
constrained  warping  and,  (b)  investigating  the  elastically  coupled  divergence 
and  flutter  of  composite  wings  with  constrained  warping,  (c)  formulating 
a  numerical  extended  beam  superelement,  which  can  capture  the  essence  of 
constrained  warping,  for  studying  geometrically  more  arbitrary  aircraft  wing 
configurations. 
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7.0  STATUS  OF  PUBLICATIONS 

The  significant  results  obtained  thus  far  in  this  research  program  are  being 
compiled  and  written  up  for  publication  in  suitable  technical  journals.  AIAA 
Paper  No.  86-1006,  a  paper  based  on  this  research  has  been  prepared  and  was 
presented  at  the  AIAA/ASME/ASCE/AHS  27th  Structures,  Structural  Dynamics  and  M 
terials  Conference  held  in  San  Antonio,  Texas,  May  19-21,  1986. 
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9.0  INTERACTIONS  WITH  AIR  FORCE  AGENCIES 

Dr.  Sinclaire  M.  Scala,  Vice  President  for  Research  and  Dr.  Gabriel  A.  Oyibo, 
the  Principal  Investigator  and  Program  Manager  for  the  research  program,  were 
invited  by  Dr.  James  Olsen,  Assistant  for  Structures  and  Dynamics  Division, 

FDL,  WPAFB  on  13  August  1985  to  brief  him  on  new  research  developments  at  FRC 
in  Aeroelasticity  and  Transonic  aerodynamics.  Dr.  Olsen  and  his  colleagues 
were  happy  with  the  briefing  which  included  the  presentation  of  some  of  the  re¬ 
sults  included  in  this  report.  They  expressed  interest  in  expanding  the  re¬ 
search  program  to  emphasize  its  application  aspects  for  the  ATF. 
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10.0  OISCOVERIES  ANO  INVENTIONS 

The  main  discoveries  of  this  research  program  thus  far  include  the  following: 
(1)  proper  modelling  of  constrained  warping  phenomenon  is  very  important  in 
aeroelastic  analysis,  (2)  this  phenomenon  is  particularly  important  for  wings 
with  coupling  and  higher  oscillation  frequencies,  (3)  there  exist  closed-form 
dynamic  solutions  for  composite  wings  with  elastic  coupling  and  warping  con¬ 
straint;  (4)  the  mechanism  for  dynamic  modal  transformations  in  composite  wings 
may  be  linked  to  some  conservative  inter-modal  energy  transfer  and 
(5)  chordwise  bending  is  more  important  for  bending  modes  than  twisting  modes 
for  or  orthotropic  wings. 
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